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1  INTRODUCTION 


1.1  Objectives 

The  technical  objectives  of  the  present  project  were: 

1.  to  obtain  experimental  data  for  heat  transfer  coefficients  for  flow  of  carbon  dioxide  (CO2, 
R-744)  in  a  micro-channel  tube,  with  an  inner  diameter  of  approximately  0.8  mm,  under 
relevant  conditions  for  compact  gas  coolers  and  evaporators  in  C02-based  air 
conditioning  systems. 

2.  to  obtain  experimental  data  for  pressure  drop  gradients  for  flow  of  carbon  dioxide  in 
micro-channel  tubes,  with  an  inner  diameter  of  approximately  0.8  mm,  under  relevant 
conditions  for  compact  gas  coolers  and  evaporators  in  C02-based  air  conditioning 
systems. 

3.  to  assess  the  consequences  of  the  measured  data  regarding  design  of  compact  heat 
exchangers  for  C02-based  air  conditioning  systems.  Of  particular  interest  is  the 
assessment  of  reversible  heat  exchangers  for  cooling/heating  operation. 

4.  to  enable  the  US  Army  to  accelerate  the  development  of  efficient  heat  exchangers  for  C02 
air  conditioning  and  heat  pump  units. 


Operation  near  the  critical  point  (31.1°C/73.8  bar)  is  expected  to  affect  heat  transfer 
characteristics  of  C02.  Earlier  measurements  indicate  effects  on  heat  transfer  from  property 
variations  (wall-bulk  variation,  buoyancy  effects)  near  pseudo-critical  points  at  super-critical 
pressures,  as  well  as  nucleate  boiling  regimes  with  improved  heat  transfer  for  evaporation  at  near- 
critical  pressures  [1].  There  is  considerable  uncertainty  about  the  validity  of  pressure  drop 
correlations  developed  for  larger  diameter  tubes  when  used  for  modelling  small-diameter  tube 
flow,  especially  during  evaporation. 


1.2  Background 

Chlorine-containing  CFC  and  HCFC  refrigerants  used  in  civil  and  military  air  conditioners  are 
now  being  phased  out  due  to  their  ozone  depleting  effect.  The  fluorocarbon  chemical  industry  is 
offering  chlorine-free  HFC  fluids  as  replacement  refrigerants,  but  these  fluids  still  have 
considerable  global  warming  impact,  and  strict  measures  will  have  to  be  introduced  to  prevent 
escape  of  refrigerant  to  the  atmosphere.  These  measures  will  include  mandatory 
systems/procedures  for  recycling  and  recovery  that  may  seriously  impact  the  logistics  in  military 
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operations.  There  is  also  some  uncertainty  about  long-term  health  and  environmental  effects  of 
these  man-made  chemicals.  Another  approach  for  developing  sustainable  CFC/HCFC-free 
technology  is  to  use  naturally  occurring  fluids  as  refrigerants  [3].  This  is  a  relatively  new  area  of 
research  -  that  is  now  gaining  interest  especially  in  Europe.  For  obvious  reasons,  these 
development  do  not  enjoy  the  support  of  the  fluorocarbon  industry. 

The  natural  refrigerant  carbon  dioxide  (CO2,  R-744-)  has  recently  been  "rediscovered"  as  a 
possible  vapour  compression  working  fluid  [4],  [5],  mainly  due  to  its  unique  combination  of 
ecological  and  personal  safety.  C02  is  a  non-toxic  and  non-combustible  natural  fluid  that  offers 
full  environmental  safety.  It  is  widely  available  in  sufficient  quantities  and  at  a  very  reasonable 
cost  all  over  the  world,  and  there  is  no  need  to  recover  or  recycle  it.  The  thermodynamic  and 
transport  properties  of  C02  are  quite  different  from  those  of  conventional  refrigerants.  As  a 
consequence,  novel  cycles  and  circuiting  concepts,  new  control  schemes  [6]  and  redesigned 
components  (compressors,  heat  exchangers)  need  to  be  developed. 

A  number  of  research  and  development  projects  on  C02-based  vapour  compression  systems  are 
now  in  progress,  focusing  on  applications  such  as  motor  vehicle  air  conditioning  systems  [4]  and 
heat  pumps  [7],  hot  water  heat  pumps  [8],  residential/commercial  air  conditioning/heat  pumps 
[9],  and  various  refrigeration  areas.  The  common  goal  for  these  projects  is  to  determine  or 
demonstrate  the  practical  potential  of  C02  systems.  The  Norwegian  University  of  Science  and 
Technology  (NTNU)  and  its  contract  research  foundation  SINTEF  has  conducted  extensive 
research  in  this  area  since  1988,  and  is  one  of  the  leading  laboratories  in  the  world  on  C02  air 
conditioning  technology. 

In  the  USA,  the  Environmental  Systems  Branch  of  the  CECOM  RD&E  Center  has  initiated 
investigations  related  to  C02  air-conditioning  units  for  the  US-Army  [10],  The  goal  is  to  develop 
smaller  and  lighter  packaged  A/C  or  HVAC  units  for  military  applications,  using  the 
environmentally  neutral  refrigerant  C02. 

The  weight  and  volume  of  heat  exchangers  built  with  MultiPort  Extruded  (MPE)  aluminium 
micro-channel  tubes  can  be  reduced  compared  to  conventional  tube/fin  heat  exchangers  with 
mechanically  expanded  round  tubes.  This  present  project  was  aimed  at  basic  research  on  heat 
transfer  coefficients  and  pressure  drop  data  inside  MPE  tubes  with  C02  as  the  refrigerant,  and  the 
resulting  impact  on  heat  exchanger  design.  The  results  of  this  project  will  be  data  that  can  be  used 
for  design  optimisation  of  MPE  tube  heat  exchangers  [11].  This  design  optimisation  will  focus  on 
heat  exchanger  performance  as  well  as  size  and  mass  reduction. 

Based  on  extensive  communication  between  NTNU/SINTEF  and  the  Environmental  Systems 
Branch  of  CECOM  RD&E  Center  (John  Manzione),  it  was  decided  to  apply  for  Seed  Project 
funding  from  ERO  to  establish  R&D  Cupertino,  and  start  the  research  necessary  for  developing  a 
new  generation  of  C02-based  Environmental  Control  Units. 
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1.3  This  report 

This  report  presents  the  results  of  the  investigations.  The  contents  are: 

•  Chapter  1 :  Basic  correlations  for  heat  transfer  and  pressure  drop 

•  Chapter  2:  Description  of  the  test  facilities  at  the  laboratory  of  SINTEF-NTNU 

•  Chapter  3:  Background  of  the  Wilson  Plot  calibration  which  is  used  to  find  a  reliable 
correlation  of  the  water-side  heat  transfer  coefficient 

•  Chapter  4:  Shows  the  results  of  the  heat  transfer  and  pressure  drop  measurements 

•  Chapter  5:  Summarises  the  main  results 

•  Chapter  6:  Gives  a  brief  assessment  of  the  implications  of  the  measurement  data  on 
evaporator  design 

The  Appendix  contains  recommended  correlations  for  heat  transfer  and  pressure  drop  of  CO2  in 
the  MPE  tube. 

1.4  Acknowledgements 

It  has  to  be  stated  that  Mr.  Andreas  Leister  has  carried  out  most  of  the  experiments  within  his 
diploma  work  at  SINTEF  Energy  Research  [1].  His  contribution  is  gratefully  acknowledged. 

The  project  was  sponsored  jointly  by  the  European  Research  Office,  London,  and  the  CECOM 
Research,  Development,  and  Engineering  Center  of  the  U.S.  Army  at  Ft.  Belvoir,  Virginia. 
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2  BASIC  HEAT  TRANSFER  AND  PRESSURE  DROP 

CORRELATIONS 


The  following  text  outlines  some  basic  equations  for  heat  transfer  and  pressure  drop  which  are 
used  in  order  to  obtain  the  experimental  results  presented  in  Chapter  5.  Here,  largely  the  same 
approach  is  used  as  in  [25]. 

The  heat  transfer  coefficient  is  defined  as 


where  q  is  the  heat  flux  according  to 


Q  is  the  heat  transfer  rate  and  .4  is  the  heat  transfer  area,  i.e.  the  tube  wall  surface,  AT  refers  to 
the  local  temperature  difference  between  the  fluid  (7)  and  the  wall  (Tw): 

AT  =  \T-TW\,  (3) 

Unfortunately  it  is  not  as  easy  as  it  seems  to  measure  this  temperature  difference  between  the 
bulk  and  the  wall  surface,  because  high  heat  transfer  coefficients  lead  to  small  temperature 
differences  to  measure.  Furthermore,  it  is  difficult  to  place  temperature  sensors  in  small  tubes, 
especially  in  a  MultiPort  Extruded  (MPE)  tube.  Therefore,  another  method  was  applied  to 
determine  the  heat  transfer  coefficient  of  the  refrigerant.  The  procedure  is  based  on  the  heat 
exchange  between  CO2  and  a  ‘secondary’  fluid  (water).  The  following  paragraphs  give  an 
introduction  to  the  heat  transfer  process  in  a  heat  exchanger. 

Figure  1  shows  a  small  cell  of  a  counter-flow  heat  exchanger  with  the  length  A l.  In  the  following, 
Al  is  assumed  to  be  small,  thus  all  fluid  properties  and  the  heat  transfer  coefficient  may  be 
considered  constant. 
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_ Length  Al _ 

1  2 

Figure  1:  A  model  of  a  small  counter-flow  heat  exchanger  cell 

The  index  ‘c’  refers  to  the  cold  fluid  and  ‘h’  to  the  hot  fluid.  The  hot  fluid  enters  the  element  at 
point  1  and  leaves  it  at  point  2  while  the  cold  fluid  enters  at  2  and  leaves  1. 

Because  of  the  first  law  of  thermodynamics  in  an  ideally  insulated  heat  exchanger  the  heat 
delivered  by  the  hot  fluid 

Q  =  Mh-{hh  v-hh2),  (4) 

must  be  equal  to  the  heat  absorbed  by  the  cold  fluid: 

Q  =  Mc-(hca-hcA),  (5) 

Now,  it  is  necessary  to  combine  the  heat  exchanged  with  the  overall  heat  transfer  coefficient: 

Q  =  k-Ah-ATlm.  (6) 

where  k  is  the  overall  heat  transfer  coefficient  based  on  the  heat  transfer  surface  on  the  hot  side 
{Ah): 


2_ 

k 


=—+K 

CCu 


cr  -A„ 


(7) 
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Ac  is  the  heat  transfer  surface  on  the  cold  side,  which  is  here  the  outer  side.  The  heat  transfer 
resistance  ( Rw )  in  the  tube  wall  is  based  on  the  inner  tube  area.  The  possible  existence  of  fouling 
layers  is  disregarded. 

As  stated  above,  a  small  heat  exchanger  element  is  considered  in  which  all  properties  and  even 
the  temperatures  are  constant.  Since  the  present  test  section  is  not  as  small,  the  temperatures 
change.  Therefore,  it  is  necessary  to  replace  the  (local)  temperature  difference  by  the  mean 
temperature  difference,  the  logarithmic  mean  temperature  difference  (A Tlm )  [14]: 


A  Tlm  = 


at,-at2 
(at, 


In 


at; 


(8) 


with 


AT;  =  Th  l  -  Tcl  and  A T2=Tha-Tc<2.  (9) 

Equation  (8)  requires  constant  properties  of  both  fluids  and  a  constant  overall  heat  transfer 
coefficient.  The  first  condition  was  not  generally  fulfilled.  In  particular  for  the  single-phase 
experiments,  the  gliding  temperature  on  the  C02-side  caused  the  specific  heat  capacity  cp  to  vary 
significantly,  especially  in  the  vicinity  of  the  (pseudo-)  critical  point.  Therefore,  the  following 
approach  was  used  instead  to  obtain  a  mean  temperature  difference  from  which  to  calculate  the 
overall  heat  transfer  coefficient:  In  the  calculation,  the  test  section  was  divided  into  60  cells  i  of 
equal  heat  exchange  Q, .  Using  the  measured  inlet  and  outlet  conditions,  as  well  as  an  assumption 
that  the  pressure  decreased  linearly  with  the  exchanged  heat  (see  [15]),  the  temperature  profile  of 
each  fluid  could  be  calculated.  Thus  the  temperature  T)  of  each  cell  was  known.  Rewriting 
equation  (6)  for  cell  i  yields 


(M\  = 


(10) 


where  At  is  the  CCVside  area  necessary  to  exchange  the  heat  in  cell  i  and  A T ,  is  the  mean 
temperature  difference  between  the  hot  and  the  cold  fluid  of  cell  i.  It  was  inherent  in  the  design  of 
the  experiment  that  the  overall  heat  transfer  coefficient  k  was  assumed  to  be  constant.  Now  it 
could  be  calculated  by  summing  equation  (10)  over  all  the  cells  N,  and  dividing  by  the  known 
total  C02-side  area,  A  a 
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(11) 


Ah  <=i 

As  shown  above,  the  overall  heat  transfer  is  given  by  the  heat  transfer  coefficients  on  the  hot  and 
the  cold  side  and  the  (known)  heat  transfer  resistance  in  the  wall.  If  one  of  the  heat  transfer 
coefficients  ( ah  or  ac)  is  known,  it  is  possible  to  determine  the  unknown  (ac  or  ah).  In  the 
present  tests,  water  flows  on  the  outer  side  of  the  tube  and  the  refrigerant  (C02)  flows  on  the 
inner  side.  In  order  to  find  an  equation  for  the  single-phase  heat  transfer  coefficient  of  the  water 
the  so-called  Wilson-Plot  method  was  used.  This  method  is  described  in  Chapter  4.  By  means  of 
this  value,  the  (re-)calculation  of  the  heat  transfer  coefficient  on  the  C02-side  can  be  easily  done 
if  the  conditions  of  both  fluids  at  the  inlet  are  known. 

The  second  part  of  the  present  work  was  to  determine  the  pressure  drop  of  the  refrigerant  in  the 
MPE  tube  along  the  heat  transferring  length.  Friction  and  acceleration/deceleration  of  the 
refrigerant  cause  this  drop.  The  design  of  the  test  section  allows  only  the  measurement  of  the  total 
pressure  drop  between  the  manifolds  and  not  the  measurement  of  the  pressure  drop  along  the  heat 
transfer  length  (Figure  2  and  see  Chapter  3).  Therefore,  it  was  necessary  to  do  some  assumptions. 
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Figure  2:  Sketch  of  the  MPE  tube  (with  water  channels)  and  refrigerant  pressure  along  the  tube 

(simplified) 

While  flowing  from  the  manifold  into  the  MPE  tube  -  consisting  of  25  parallel  ports  -  a  sudden 
contraction  occurs  to  the  refrigerant.  This  leads  to  a  pressure  drop  at  the  inlet.  At  the  outlet  of  the 
MPE  tube  the  opposite  effect  occurs.  An  expansion  leads  to  a  pressure  rise.  Calculation  models 
for  these  pressure  changes  are  described  in  [17]  and  [19]  ([15]): 

In  case  of  a  single-phase  flow  the  contraction  pressure  drop  at  the  inlet  is  given  by 


(12) 
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The  pressure  rise  due  to  expansion  at  the  outlet  can  be  calculated  according  to 


A Po 


(13) 


Where  m  is  the  mass  flux  of  the  fluid  in  the  ports  and  a  is  the  ratio  of  the  total  port  area  ( APytot ) 
and  the  front  area  (Ajj)  of  the  MPE  tube.  For  the  present  test  tube  it  is 


G  = 


0.28 


(14) 


The  entrance  and  exit  pressure  loss  coefficients  (Kc  and  Ke)  are  also  a  function  of  the  contraction 
and  expansion  geometries  and  in  some  cases  of  the  Reynolds  number.  Kays  and  London 
determined  them.  For  turbulent  flow  with  an  approximate  Reynolds  number  of  10,000  and  the 
geometry  of  the  used  tube  these  factors  are 


Kc  =  0.4  and  Ke  =  0.5 .  (15) 

In  case  of  a  two-phase  flow,  the  models  for  the  effects  at  the  inlet  and  the  outlet  are  more 
complicated  ([18], [19]).  The  inlet  pressure  drop  caused  by  the  sudden  contraction  can  be  written 
as 


A PKC  = 


^  •  \2 

m 


VC‘J 


(1-0- 


(l  +  C,)- 

(  v3  2  /,  \3  2 

X  .  (I"*)  -V/ 

l  v  J  ; 

2- 

i*-vs+(l-x)-v,J 

f  ..2 


C  - 


M) 


(16) 


and  the  outlet  pressure  rise: 


A p0,e  =m-cr-{l-(j)-vl  - 


(l-*)2 

(M 


+ 


V  '  J 


(17) 


As  can  be  seen  from  the  equations  (16)  and  (17),  there  are  two  new  parameters,  the  coefficient  of 
contraction  (Cc)  and  the  void  fraction  (e): 


Cc  depends  on  the  contraction  ratio  already  defined  by  equation  (14).  For  the  present  installation, 
the  coefficient  becomes 
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Cc  =  0.6 


(18) 


The  void  fraction  can  be  calculated  as 


(pA 

2  " 
3 

£  = 

1  + - 

X 

[aJ 

The  refrigerant  also  flows  through  two  short  sections  at  the  ends  of  the  tube  where  no  heat  is 
transferred.  The  frictional  pressure  drop  occurring  in  these  adiabatic  parts  of  the  MPE  tube  may 
be  calculated  by  means  of  the  Blasius  friction  factor 


/  = 


0.3164 


Re025  ' 

and  the  general  form  for  the  pressure  drop: 


A Pf  =f— 7 


A  /  m 


d  2  p 


(20) 


(21) 


In  case  of  two-phase  flow,  the  adiabatic  pressure  drop  may  also  be  calculated  using  the  Blasius 
friction  factor  but  using  an  equivalent  mass  flux  for  the  liquid  phase  ( meq )  instead  of  the  total 
mass  flux  ( m )  [20]: 


meq=m-\ 


(l-x)+. 


le$ 

KP*J 


(22) 


With  these  equations  the  ‘adiabatic’  pressure  drop  at  the  inlet  (  A pi  f )  and  at  the  outlet  ( A po  f ) 
can  be  calculated. 

Finally,  the  pressure  drop  occurring  along  the  heat  transferring  length  of  the  test  section  (A pts ) 
can  be  calculated  according  to 

A Pts  =  A Pmeas  ~  A Pi,c  ~  A?,,/  "  AP0,/  +  A Po,e  •  (23) 

Note  that  because  of  the  horizontal  installation  of  the  MPE  tube,  there  is  no  influence  of  gravity. 
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3  TEST  FACILITIES 


This  chapter  gives  an  overview  of  the  test  rig  for  measuring  heat  transfer  and  pressure  drop  data 
of  C02  in  MPE  tubes  at  the  laboratory  of  Sintef  -  NTNU  in  Trondheim,  Norway. 

3.1  Overview 

The  test  rig  consisted  of  three  main  parts: 

•  the  refrigerant  loop, 

•  the  test  section,  and 

•  the  “coolant”  (water)  loop  1 

Figure  3  shows  the  flow  circuit  of  the  test  rig  that  allowed  the  analysis  of  the  heat  rejection 
process  of  C02  at  super-critical  pressures  (single-phase  cooling,  ‘gas  cooling’)  and  the  heat 
absorption  process  at  sub-critical  pressures  (evaporation)  in  MPE  tubes  (Figure  7,  Figure  8).  The 
refrigerant  flow  through  the  system  depended  on  the  process  to  be  analyzed.  Figure  4  a)  shows 
the  refrigerant  flow  at  cooling  tests  and  Figure  4  b)  at  evaporation  tests.  The  following  system 
description  is  based  on  the  cooling  tests.  Afterwards  the  main  differences  for  the  evaporation  tests 
are  pointed  out  briefly. 

Cooling  Tests 

For  testing  super-critical  C02,  the  desired  pressure  in  the  test  rig  was  set  with  an  electrically 
heated  C02  cylinder.  This  cylinder  was  connected  to  the  system  at  (a).  Refrigerant  was  circulated 
through  the  system  by  means  of  a  gear  pump  (1).  It  entered  the  pre-heat  section  (4)  where  heat 
was  added  electrically,  thus  the  refrigerant  state  could  be  set  to  the  desired  test  conditions,  i.e. 
heated  up  to  the  desired  temperature.  The  refrigerant  then  entered  the  test  section  (9)  where  heat 
exchange  between  the  refrigerant  and  the  coolant  took  place.  After  exiting  the  test  section,  the 
refrigerant  entered  the  refrigerant  condenser  (6),  where  was  removed. 

Figure  5  shows  the  measurement  points  of  the  refrigerant  loop,  which  included  a  Coriolis-type 
mass  flow  meter  ( Mref  ),pressure  and  temperature  measurements  at  the  preheat  section  inlet 
(trefips,  Prefips),  temperature  between  the  upper  and  lower  preheat  section  (fref.mps),  temperature  and 
the  pressure  at  the  test  section  inlet  (trefits,  pref,usX  temperature  at  the  test  section  outlet  (tref0ts),  and 
pressure  drop  across  the  test  section  (A prej).  The  cooling  fluid  temperature  (tCOn,0)  was  also 
measured. 

1  Although  water  is  used  as  heat  source  and  as  heat  sink  the  expression  coolant  loop  is  used  here. 
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Refrigerant  loop: 

1  Refrigerant  pump 

2  Refrigerant  filter 

3  Refrigerant  mass  flow  meter 

4  Pre-heat  section 

5  Safety  valve  (102  bar) 

6  Refrigerant  condenser 


7  Refrigerant  receiver  with  sight  glass 

8  ‘Sub-cooler’ 

Test  section: 

9  Test  section 
Coolant  loop: 

10  Coolant  cooler 

1 1  Coolant  tank  with  heat  exchanger 

12  Water  heater  (Type  Haake) 


13  Main  coolant  pump 

14  Booster  coolant  pump 

15  Coolant  filter 

16  Coolant  mass  flow  meter 
Supply  lines: 

a  Refrigerant  charge 
b  Coolant  charge 
c  City  water  inlet 
d  City  water  outlet 


Figure  3:  Flow  chart  of  the  MPE  tube  test  rig 


a)  refrigerant  flow  at  gas  cooling  tests  b)  refrigerant  flow  at  evaporation  tests 

Figure  4:  Flow  of  the  refrigerant  through  the  system  at  cooling  and  evaporation  tests 


Figure  5:  Refrigerant  loop  with  measurement  points  (component  no.  see  Figure  3) 


The  coolant  loop  with  the  measurement  points  is  shown  in  Figure  6.  As  can  be  seen,  after  exiting 
the  test  section  (9)  the  heated  coolant  flowed  through  a  heat  exchanger  (10)  where  it  was  cooled 
down.  In  the  following  coolant  tank  (11)  the  coolant  could  be  heated  to  a  defined  temperature 
with  a  water  heater  of  the  type  Haake  (12),  thus  the  desired  coolant  temperature  could  be  set. 
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Two  centrifugal  pumps  (13  and  14)  circulated  the  coolant  through  the  loop  and  by  means  of 
varying  the  opening  of  the  bypasses  the  flow  rate  could  be  varied. 


Figure  6:  Coolant  loop  with  measurement  points  (component  no.  see  Figure  3) 


The  mass  flow  rate  (Mcool)  was  measured  with  a  Coriolis-type  mass  flow  meter  (16). 
Temperatures  and  pressures  at  the  test  section  inlet  and  outlet  (tcooijts;  tcooi,ots',  Pcooiitsl  Pcooiots) 
were  also  measured.  Close  to  the  mass  flow  meter,  the  temperature  of  the  coolant  {tcooi,m)  was 
measured.  Another  measured  temperature  was  the  coolant  temperature  at  the  cooler  inlet  ( tcooier,i) ■ 


Evaporation  Tests 

The  main  adjustment  necessary  to  run  the  evaporation  tests  was  the  change  of  the  cooling  fluid 
used  in  the  CO2  condenser  (6),  because  temperatures  below  0°C  were  necessary  (down  to  about 
-15°C).  Ethanol  cooled  by  a  Haake  COOLER  was  used  for  evaporation  tests  instead  of  city  water. 

Refrigerant  flow  through  the  system  in  evaporation  tests  was  also  different.  It  is  shown  in  Figure 
4  b).  To  obtain  the  desired  saturation  pressure  of  the  refrigerant,  the  refrigerant  leaving  the 
condenser  flowed  through  the  vessel  (7)  before  it  reached  the  pump.  The  refrigerant  charge  had 
to  be  chosen  in  such  a  way  that  during  the  tests,  a  liquid  level  could  be  seen  in  this  vessel,  i.e.  the 
refrigerant  left  the  condenser  in  saturated  liquid  state.  To  avoid  cavitation  in  the  refrigerant  pump 
located  after  the  vessel,  the  refrigerant  had  to  be  cooled  down  (sub-cooling  of  about  3  K).  For 
this  purpose  a  small  heat  exchanger  (8)  was  installed. 
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3.2  Components 


3.2.1  Test  Section 

The  MPE  tube  to  be  investigated  was  a  540  mm  long  tube  soldered  in  manifolds  on  both  ends. 
The  cross-section  of  the  MPE  tube  can  be  seen  in  Figure  7.  The  tube  was  made  of  a  zinc-coated 
aluminium  tube  with  25  round  ports  each  with  a  diameter  of  0.787  mm  (manufactured  by  Hydro 
Aluminium  Adrian). 


Figure  7:  Cross-section  of  the  multiport  extruded  (MPE)  tube 


Figure  8:  Photograph  of  the  MPE  test  tube  with  manifold  lying  in  the  (opened)  Teflon  jacket 

The  test  tube  was  placed  inside  a  PFTE  (Teflon)  jacket,  where  the  coolant  channels  were  milled 
out.  These  channels  were  shaped  like  an  edgy  double  helix  around  the  MPE  tube.  A  photograph 
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of  the  MPE  tube  laying  in  the  (opened)  PFTE  jacket  is  shown  in  Figure  8.  The  outside  of  the 
jacket  was  insulated  with  25  mm  thick  plates  of  polystyrene  (Styrofoam). 

3.2.2  Refrigerant  Pump 

The  refrigerant  pump  (1)  was  a  gear  pump  of  the  type  MicroPump  219.  By  means  of  controlling 
the  voltage  (0  to  24  V)  it  was  possible  to  vary  the  number  of  revolutions  infinitely  between  0  and 
3,000  rpm.  The  maximum  capacity  of  the  pump  was  2.5  1/min. 

3.2.3  Coolant  Pumps 

The  main  pump  (13)  is  a  centrifugal  pump  of  the  type  GRUNDFOS  CHI  4-60  with  a  capacity  of 
10 1/min  at  6  bar.  It  was  driven  by  a  220  V  AC  motor. 

During  the  Wilson  Plot  calibration  and  the  experiments  it  as  necessary  to  vary  the  coolant 
(water)  flow  rate  in  a  wide  range.  In  the  high  end  of  the  range  a  large  pressure  drop  had  to  be 
overcome.  Therefore  a  ‘booster’  pump  was  needed  (14).  This  pump  was  a  centrifugal  pump  of 
the  type  Grundfos  CH  2-50  with  a  capacity  of  10 1/min  at  4  bar.  It  was  connected  in  series  to 
the  main  pump.  This  pump  was  only  used  when  high  coolant  flow  rates  were  required.  The 
maximum  coolant  flow  rate  was  about  12 1/min  when  both  pumps  were  used. 

3.2.4  Bypasses 

To  regulate  the  coolant  flow  rate  and  to  discharge  the  heat  delivered  by  the  coolant  pumps 
several  bypasses  were  placed  in  the  coolant  loop.  The  main  bypass  was  from  the  outlet  of  the 
booster  pump  (14)  to  the  inlet  of  the  coolant  cooler  (10).  This  was  used  to  ‘discharge’  the  heat 
delivered  by  the  pumps  to  the  coolant.  While  the  pumps  were  running,  the  bypass  could  be 
controlled  with  a  valve  such  that  the  temperature  did  not  increase.  To  control  the  flow  rate 
through  the  coolant  loop,  there  was  a  valve  at  the  outlet  of  the  booster  pump  (14).  To  adjust  the 
mass  flow  rate  more  precisely,  two  bypasses  were  placed  around  this  valve.  There  was  also  a 
bypass  around  the  coolant  filter  (15),  because  a  huge  pressure  drop  occurred  over  the  filter  in 
case  of  high  flow  rates.  When  a  high  flow  rate  was  needed,  this  bypass  could  be  opened. 

3.3  Measurement  Equipment 

This  section  describes  the  instrumentation  and  measurement  components  used.  During  the  tests  it 
was  necessary  to  measure  temperatures,  pressures,  mass  flow  rates,  and  electrical  power.  All 
these  data  were  logged  and  processed. 

3.3.1  Temperature 

Two  different  methods  were  used  to  measure  the  temperatures.  Some  temperatures  were 
measured  with  thermocouples,  and  the  more  important  temperatures  (in  the  test  section)  were 
measured  with  Pt  100  sensors. 

The  thermocouples  were  of  the  type  T  ,  Copper  (Cu)  /  Constantan  (CuNi).  The  reference 
temperature  was  given  by  crushed  ice.  Pt  100  elements  measured  the  temperatures  required  to 
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determine  the  heat  transfer  rate  in  the  test  section,  i.e.  the  temperatures  at  the  inlet  and  the  outlet 
of  the  test  section  on  the  refrigerant  side  (trefus->  Uefots )  and  on  the  coolant  side  (tcooi,  its -i  tcool,ots)- 


3.3.2  Pressure 

The  pressures  prefiPs  and  prepts  were  each  measured  with  a  pressure  transducer  of  the  type 
Honeywell  STA3000  (gauge  =  pressure  over  atmospheric  pressure.  The  atmospheric  pressure 
of  the  laboratory  was  added  to  obtain  absolute  pressure).  The  refrigerant  pressure  difference 
between  test  section  inlet  at  outlet  (A pref)  was  measured  with  a  HONEYWELL  STD924  transducer. 
Manometers  were  also  connected  to  the  coolant  inlet  and  outlet  to  determine  the  pressure  drop  on 
the  coolant  side  along  the  test  section. 

3.3.3  Mass  Flow 

Refrigerant  mass  flow  rate  was  measured  with  a  Coriolis-type  Danfoss  Mass  2100  sensor 
connected  to  a  Danfoss  Mass  3000  signal  converter. 

Coolant  mass  flow  rate  was  also  measured  with  a  Coriolis-type  mass  flow  meter,  but  the  sensor 
was  a  Rheonik  RHM  04  GNT  connected  to  a  Rheondc  RHE  08  signal  converter. 

3.3.4  Electrical  Power 

The  electrical  power  for  the  pre-heat  section  was  regulated  by  means  of  two  variable  resistances 
and  displayed  on  two  ELCONTROL  VIPD  watt  transducers.  The  logged  signals,  however,  were 
supplied  by  two  L-Unit  LWT-24-A1-H  watt  transducers. 

3.3.5  Data  Logging  and  Processing 

Data  from  thermocouples,  pressure  transducers,  mass  flow  meters  and  watt  transducers  were 
processed  and  logged  by  a  Keithley  7001  Switch  System  and  a  Keithley  2010  Multimeter. 
The  switcher  changed  between  the  different  input  channels,  while  the  multimeter  measured  the 
voltage,  digitised  it,  and  sent  it  to  a  PC. 

The  Keithley  data  logger  was  not  capable  of  measuring  with  Pt  100-elements.  Therefore  an 
ASA  F250  MK II  Precision  Thermometer  was  used  combined  with  an  A2A  SB250  Multichannel 
switch  box. 

The  PC  processed  the  data  by  a  computer  program  developed  at  Sintef  Energy  Research.  It 
converted  the  data  to  values  with  physical  units.  The  data  could  also  be  shown  as  graphs  on  the 
computer  screen.  At  the  same  time  the  program  wrote  them  to  a  file  which  could  be  processed  by 
other  programs. 

After  the  measurement,  these  data  files  could  be  analysed  in  MS  Excel97  worksheets.  The 
thermodynamic  and  transport  properties  of  the  refrigerant  and  of  the  coolant  were  supplied  by 
the  libraries  xlco21ib.dll  and  h2oprop_p.xls,  respectively,  which  were  also  developed  at 
Sintef  Energy  Research. 
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xlco21ib.dll  is  an  implementation  of  the  IUPAC  thermodynamic  equations  [42]  with 
improvements  [43],  as  well  as  transport  properties  [44]. 
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4  WILSON  PLOT  CALIBRATION 


As  mentioned  in  Chapter  2,  it  is  necessary  to  have  knowledge  of  the  heat  transfer  coefficient  on 
the  ‘coolant’  (water)  side  (ac)  to  determine  the  heat  transfer  coefficient  of  the  refrigerant  (ah) 
by  means  of  the  heat  transfer  rate  (Q)  and  the  overall  heat  transfer  coefficient  (k)  in  a  heat 
exchanger.  On  beforehand,  ac  was  unknown.  Common  correlations  for  turbulent  single-phase 
flow  could  not  be  used,  due  to  the  geometry  of  the  test  section  and  the  desired  accuracy  of  the 
data.  A  correlation  for  the  water-side  heat  transfer  coefficient  ac  was  found  using  a  modified 
Wilson  plot  method 

A  short  description  of  the  principles  can  be  found  in  Section  4.1.  The  results  are  presented  in 
Section  4.2. 


4.1  Principles  of  the  Wilson  Plot  Method 

The  Wilson  plot  method  was  introduced  in  1915  by  E.E.  Wilson  as  a  tool  for  designing  shell- 
and-tube  condensers  and  feed  water  heaters  [23].  In  this  work,  the  Wilson  plot  method  is 
modified  as  described  in  [25]. 

The  Dittus-Boelter-type  equation  is  a  well-known  expression  for  the  single-phase  heat  transfer 
coefficient  in  a  turbulent  flow: 


a  c=C-Rem-Prn-—  (24) 

dhycl 

This  equation  is  used  to  describe  the  heat  transfer  coefficient  on  the  water-side.  For  this  purpose, 
the  Reynolds  number  exponent  (m)  and  the  constant  (C)  are  assumed  to  be  unknown,  but  they 
can  be  found  with  the  Wilson  Plot  technique.  The  Prandtl  number  exponent  (n)  is  commonly 
taken  to  be  as  n  =  0.4  for  heating  and  n  =  1/3  for  cooling  [14]. 

Combining  the  equation  for  the  overall  heat  transfer  coefficient  (7)  with  equation  (24)  yields 
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This  is  an  equation  on  the  form 


(25) 


7,  =A-Xt  +B. 


(26) 
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where  the  slope  A  equals 


A=- 


the  independent  variable  is 
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and  the  intersection  with  they-axis  (B)  equals 


(27) 


(28) 
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The  dependent  variable  ( Yj )  is  defined  as 


(29) 


r,= 


I'* 


(30) 


In  Equation  (30),  the  overall  heat  transfer  coefficient  (k)  is  known  from  the  experiment,  since  the 
transferred  heat,  the  heat  transferring  area  and  the  mean  temperature  difference  is  given  by  the 
measurements  and  the  geometry.  The  wall  resistance  Rw  may  be  calculated  (or  neglected)  and  Ac, 
Ah,  and  dhyd  are  geometry  data.  Thus  Y;  is  known. 

The  Reynolds  number  (Re)  and  the  Prandtl  number  (Pr)  as  well  as  the  thermal  conductivity  (/l) 
are  known  from  the  coolant  properties.  Since  m  is  not  known  in  Equation  (28),  Xj  is  unknown. 
During  the  experiments,  the  independent  variable  Xj  can  be  varied  by  changing  the  mass  flow  of 
the  coolant  (Mc).  This  corresponds  to  a  variation  of  the  Reynolds  number.  In  all  experiments, 
the  coefficient  B  has  to  be  constant,  i.e.  ah  has  to  be  constant  during  the  plot.  Therefore,  the 
following  conditions  are  required: 

•  The  mass  flow  Mh  had  to  be  constant. 

•  The  hot  side  mean  temperature  th  and  the  inlet  and  outlet  temperatures  (trefus  and  tref0ts )  had 
to  be  constant  (constant  heat  transfer  rate  Q). 

•  In  order  to  keep  Q  constant  while  varying  Mc ,  the  coolant  inlet  and  outlet  temperature,  i.e. 

A Tc ,  and  the  mean  temperature  (tc)  had  to  be  varied. 


With  the  experimental  results,  a  plot  of  Equation  (25)  can  be  made.  Figure  9  shows  such  a  plot. 
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The  line  PjP.,  is  a  linear  curve  fit  to  the  measured  data.  If  it  is  assumed  that  the  heat  transfer 
coefficient  on  the  coolant  side  becomes  infinitely  high  for  an  infinite  coolant  mass  flow  ( M c ), 
the  line  may  be  extrapolated  to  the  y- axis.  The  intersection  is  indicated  with  point  P3.  This  point 
corresponds  to  an  infinite  Reynolds  number  ( Rec )  and  therefore  to  an  infinite  coolant  mass  flow 
rate,  thus  X\  becomes  zero  (see  Equation(28)).  The  length  OP3  then  represents  the  thermal 
resistance  on  the  inner  side  of  the  tube  ( B ). 

Solution  Algorithm 

In  Equation  (25)  there  are  three  unknowns:  <%,  C,  and  m.  Unfortunately  it  is  not  possible  to  solve 
a  linear  plot  with  more  than  two  unknowns.  Therefore  a  solution  algorithm  must  be  found  (curve 
fit).  As  can  be  seen  in  Figure  11,  this  algorithm  consists  mainly  of  two  steps.  First  a  value  of  m 
has  to  be  guessed,  and  then  the  plot  can  be  made.  With  a  second  plot,  the  so-called  logarithmic 
Wilson  Plot  (Figure  10),  the  ‘goal  seek’  for  m  can  be  checked: 


Equation  (25)  can  be  written  as 
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(31) 


or 
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Taking  the  logarithm  of  (32)  results  in 


(32) 
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In  y2  =  In 


or 


In  y2  =  -m  •  In  (i?e)- In  C 
This  equals  also  to  a  linear  equation  of  the  form 


Y2=D-X2+E. 

where  the  slope  D  is  defined  with 
D  =  -m 

and  the  independent  variable  X2  is  defined  as 
X7=\a(Re) 

as  well  as  the  intersection  with  the  jy- axis  (E)  equals 
E  -  -In  C . 

By  re-substituting  >>2  we  get  the  dependent  variable  Y2  as 


(33) 


(34) 


(35) 


(36) 


(37) 


(38) 


If  the  m  from  equation  (36)  and  the  initially  guessed  m  are  the  same,  then  the  solution  for  m  is 
found.  If  not,  a  more  qualified  Reynolds  exponent  m  has  to  be  guessed.  In  Microsoft  Excel97 
this  can  easily  be  done  by  using  the  SOLVER  function.  With  this  algorithm,  one  can  find  the 
constant  (C)  and  the  Reynolds  number  exponent  (m)  needed  for  the  Dittus-Boelter-type  equation 
(24).  See  Figure  11. 

It  was  assumed  that  this  equation  described  the  heat  transfer  coefficient  on  the  outer  side  of  the 
tube  as  long  as  the  flow  was  turbulent. 
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4.2  Results  of  the  Wilson  Plot 


4.2.1  Water  Heating 

During  this  Wilson  Plot,  the  refrigerant  was  supercritical  with  a  pressure  of  pco2  =  (91  ±  0.9)  bar 
and  an  inlet  temperature  of  tco2,its=  (30  ±  0.2)°C.  The  heat  flux  was  kept  constant  at  q- 
(20  ±  0.65)  kW/m2.  Since  the  accuracy  of  the  calibration  is  improved  with  a  higher  heat  transfer 
coefficient  on  the  refrigerant-side  ( ah ),  a  mass  flux  of  mC02  =  (2000  ±  42)  kg/(m2-s)  was  chosen. 
These  data  corresponded  to  a  temperature  drop  of  ATco2  =  (7.7  ±  0.2)  K  in  the  test  section. 

The  mean  temperature  in  the  test  section,  tC02-  (25.7  ±  0.4)°C,  was  close  to  the  room 
temperature  in  the  laboratory;  therefore  the  heat  loss  to  the  ambient  was  minimal. 

The  measurements  were  made  when  the  desired  values  were  reached  and  the  rig  was  operating  in 
steady-state  conditions.  All  data  were  logged  and  written  to  a  file  by  the  PC.  For  each  measured 
point,  about  20  readings  were  taken  with  a  time  interval  of  60  seconds.  The  calibration  was  made 
with  a  total  of  29  different  water  mass  flow  rates  ( MH20  =  2  to  9.8  kg/min). 

By  means  of  the  Wilson  plot  technique,  the  graphs  shown  in  Figure  12  and  Figure  13  were 
found.  The  plot  gave 

C  =  0.202  and  m  =  0.607.  (40) 

and  the  ‘calibrated’  equation  for  the  water  side  heat  transfer  coefficient  (ac)  during  heating 
mode  had  the  form 

<xc,to  =  0.202  •  Re0607  Pr04  ■\/ dhyd .  (41) 


rrto2  =  2000  kg/(m2-s);  qco2  =  20,000  W/m2 ;  tco2  =  25  °C 


Xi  =  l/(Rem-Pr”-?Jdh) 

Figure  12:  Wilson  plot  for  water  heating  (C  =  0.202  and  m  =  0.607) 
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mco2  =  2000  kg/(m2-s);  qco2  =  20,000  W/m" ;  tco2  =  25  °C 


As  described  in  Section  4.1,  C  is  the  inverse  of  the  slope  of  the  fitted  curve  (straight  line)  in 
Figure  12  and  m  is  the  negative  of  the  slope  of  the  fitted  curve  (straight  line)  of  Figure  13. 

4.2.2  Water  Cooling 

It  was  decided  to  perform  a  second  Wilson  plot  calibration  for  the  evaporation  tests  when  the 
water  was  cooled  by  the  C02.  For  this  calibration,  the  test  conditions  on  the  inside  of  the  tube 
(refrigerant  side)  again  had  to  be  constant.  The  C02  pressure  was  set  to  pC02  =  (91  ±  0.7)  bar  and 
the  test  section  inlet  temperature  was  constant  at  tco2,us  —  (4.7  ±  0.1)°C.  With  a  heat  flux  of  q  = 
(20  ±  0.38)  kW/m2  and  a  specific  mass  flux  of  mc02=  (2000  ±  14)  kg/(m2-s),  the  resulting 
temperature  rise  and  the  mean  refrigerant  temperature  was  A Tco2=  (10.6  ±  0.2)  K  and  tC02- 
(11.3  ±  0.5)°C,  respectively.  Experiments  at  a  total  of  20  different  water  mass  flow  rates  were 
carried  out  ( MH2Q  =  1.7  to  9.4  kg/min). 

The  graphs  of  the  Wilson  Plot  are  shown  in  Figure  14  and  Figure  15,  and  the  resulting  Wilson 
Plot  constants  were 

C  =  0.162  and  m  =  0.653  (42) 

The  resulting  equation  for  the  water  heat  transfer  coefficient  in  cooling  mode  was 

CW  =  0.162 ■Re°‘,s  ■ Pr «3  , X/d (43) 
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Xi  =  l/(Rem*Prn*A/dh) 


Figure  14:  Wilson  Plot  for  water  cooling  (C  =  0.162  and  m  =  0.653) 


X,  =  l/(Rem-Prn-X/4) 

Figure  15:  Logarithmic  Wilson  Plot  for  water  cooling  (C  =  0.162  and  m  —  0.653) 


4.2.3  Comparison  of  the  Wilson  Plots 

The  heating  and  cooling  Wilson  plots  did  not  show  the  same  result,  i.e.  two  different  equations 
were  yielded.  Figure  16  shows  the  calculated  heat  transfer  coefficient  according  to  the  equations 
(41)  and  (43)  at  a  temperature  of  20°C  depending  on  the  mass  flow  rate. 
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water  mass  flow  Mh2o  [kg/m  in] 

Figure  16:  Heat  transfer  coefficient  of  water 


As  can  be  seen,  the  heat  transfer  coefficient  is  slightly  higher  for  water  cooling  (evaporation  tests 
see  Chapter  5.1)  than  for  water  heating  (gas  cooling  tests  see  Chapter  5.2).  However,  most  of  the 
experiments  were  carried  out  at  water  flow  rates  between  MH20  =  3  to  5  kg/min,  where  the 
deviation  between  the  two  curves  is  quite  small. 
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5  RESULTS 


This  chapter  describes  the  experimental  results  for  the  heat  transfer  coefficient  and  pressure  drop 
at 

•  Super-critical  (cooling),  and 

•  Sub-critical  conditions  (evaporation). 

After  a  brief  description  of  the  test  conditions,  the  measured  values  are  compared  with  ‘common’ 
correlations  from  the  literature.  Since  the  heat  exchanger  used  has  a  certain  length,  mean  values 
rather  than  local  values  were  determined.  The  exact  procedure  is  described  in  Chapter  2. 


For  the  comparisons  to  calculation  models,  the  ‘local’  values  with  respect  to  the  mean  properties 
along  the  test  section  were  used.  The  heat  transfer  and  pressure  drop  data  are  compared  in  a 
quantitative  and  qualitative  way.  The  quantitative  comparisons  were  done  by  means  of  the 
average  deviation  ((Javg)  and  the  mean  deviation  (crmean): 


a 


avg 


a 


meas 


(44) 


a 


mean 


® meas 


(45) 


where  n  is  the  total  number  of  measurements,  owe  is  the  calculated  local  and  (Was  is  the 
measured  mean  heat  transfer  coefficient.  For  the  pressure  drop,  the  analogous  deviations  were 
used. 


5.1  Supercritical  Conditions  (Cooling) 

5.1.1  General 

Figure  17  shows  a  schematic  diagram  of  a  cooling  ‘test  point’.  (See  Figure  3  and  Figure  4  for  the 
flow  chart  of  the  test  rig.) 
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Figure  17:  Schematic  temperature-entropy  diagram  of  a  gas  cooling  process 

Refrigerant  entered  the  preheat  section  at  a  super-critical  pressure  and  a  sub-critical  temperature 
(1).  Then  heat  was  added  electrically  in  the  preheat  section  before  the  CO2  entered  the  test 
section  (2).  The  heat  transfer  to  the  water  cicuit  took  place  in  the  test  section.  The  cooled  CO2 
exited  the  test  section  (3)  and  was  cooled  down  in  the  CO2  cooler  before  it  reached  the 
refrigerant  pump  inlet  (1). 

In  order  to  allow  an  analysis  of  the  heat  transfer  and  pressure  drop  characteristic,  a  total  of  six 
test  series  were  carried  out.  The  tests  were  done  at  different  mass  flow  rates  (mass  fluxes,  m ), 
pressures  (p)  and  heat  fluxes  (9 ). Table  1  summarises  the  test  conditions. 


Table  1:  Overview  of  the  cooling  test  conditions 


(kW/m2)  J| 

p  (bar) 

m  (kG/Viif-  s))""- 

BBS 

■  >  ■* ■;  ' V 

-'v.  ’I  :  r:' 

1  AT  ‘  ' 

1U1  . 

H 

- 

Gas  cooling  No.  6 

- 

Gas  cooling  No.  4 

Gas  cooling  No.  3 

Gas  cooling  No.  5 

- 

Gas  cooling  No.  2 

- 

- 

Gas  cooling  No.  1 

- 

Table  2  shows  the  test  conditions  and  the  experimental  tolerances  of  the  points.  As  can  be  seen, 
the  mass  flux  was  kept  within  a  range  of  about  ±2%  and  the  pressure  showed  a  tolerance  of  ±1%. 
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Table  2:  Test  conditions  and  experimental  tolerances 


Test 

q 

(kW/m2) 

m 

(kg/(m2-s)) 

P 

(bar) 

Gas  cooling  No.  1 

1Q+1.0 

600  ±11 

91  ±0.96 

Gas  cooling  No.  2 

20^7 

1200  ±16 

91  ±0.87 

Gas  cooling  No.  3 

20 +1'5 

Z'J-0.6 

900  ±13 

91  ±  0.59 

Gas  cooling  No.  4 

20!^ 

900  ±19 

81  ±0.30 

Gas  cooling  No  5 

20+0.6 

zu-0.6 

900  ±6 

101  ±0.36 

Gas  cooling  No.  6 

20 +1'4 
zu-1.0 

600  ±12 

91  ±0.44 

Only  the  heat  flux  showed  a  wider  range,  but  this  is  of  secondary  importance,  as  also  shown  by 
the  results  (see  below). 

Each  of  the  series  consisted  of  about  five  points  at  different  mean  temperatures  (F  ).  One  point 
was  always  close  to  the  pseudo-critical  temperature  (maximum  cp)  and  at  least  two  were  above 
and  two  below  this  important  temperature.  The  pseudo-critical  temperatures  and  the 
corresponding  maximum  specific  isobaric  heat  capacities  at  the  relevant  pressures  are  given  in 
Table  3. 

The  choice  of  test  points  below  and  above  the  pseudo-critical  temperature  was  done  such  that  the 
enthalpy  differences  between  the  points  were  approximately  constant.  However,  some  series 
were  influenced  by  constraints  given  by  the  test  rig.  In  order  to  set  the  inlet  temperature  to  the 
desired  value,  the  refrigerant  was  heated  electrically  in  the  preheat  section,  but  the  capacity  was 
limited  at  Pei  =  3  kW.  This  effected  for  example  the  test  ‘Gas  cooling  No.  2’  with  m  = 
1200  kg/(m2s),  where  the  maximum  capacity  of  the  preheat  section  was  reached  at  a  rather  low 
inlet  temperature  of  about  48°C. 


Table  3:  Pseudo-critical  temperature  and  maximum  specific  heat  capacity 


-  P 
(bar) 

tpseudo 

(°C) 

Cp 

(kJ/(kg-K)) 

81 

35.2 

29.4 

91 

40.6 

12.0 

101 

45.6 

7.8 

The  temperature  drop  along  the  test  section  is  given  by  the  mass  flow  rate,  i.e.  by  the  mass  flux 
and  the  cross  section,  and  the  heat  flux  (q)  that  is  based  on  the  perimeter  of  the  ports. 


5.1.2  Heat  transfer  at  super-critical  pressures 

In  this  section,  the  measured  single-phase  heat  transfer  coefficients  are  compared  with 
calculation  models.  The  following  models  are  used: 
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>  Different  forms  of  Gnielinski’s  correlation: 

•  simple  Gnielinski:  Gnielinski’s  correlation  for  the  Nusselt  number  using  the  Haaland 

friction  factor  [21]  but  neglecting  the  influence  of  the  wall 
temperature  (Roughness  e  =  1.0  •  10'6  m). 

•  Gnielinski:  Gnielinski’s  correlation  using  the  Haaland  friction  factor  and 

considering  the  influence  of  the  wall  temperature, 

•  Gnielinski  (VDI):  Gnielinski’s  correlation  using  the  Filonenko  friction  factor  (according 

to  [28]). 

>  The  ‘special’  models  developed  for  heat  transfer  at  super-critical  pressures: 

•  Polyakov  [29] 

•  Ghajar  &  Asadi  [30]2 

>  And  the  Dittus-Boelter  correlation  [14]  with  different  Prandtl  exponents  («)  (compare 
equation  (24)  with  C  =  0.023,  m  =  0.8): 

•  n=  1/3:  which  is  the  ‘correct’  value  for  a  flow  under  cooling  mode 

•  n  =  0.4:  which  is  normally  used  for  heating  mode. 

The  comparison  of  the  experimental  data  with  these  models  showed  a  quite  good 
correspondence.  Table  4  contains  the  average  and  the  mean  deviations  of  the  32  points  measured 
and  Figure  18  to  Figure  20  visualise  the  good  correspondence  graphically. 


Table  4:  Deviation  of  the  measured  heat  transfer  to  calculation  models 


•  Jglp:  ^ 

References 

CTavg  [%] 

Omean  [%*] 

see  above 

-3 

4 

|H  %  ■  Gnielinski  ■  1 

see  above 

-1 

4 

[28] 

-7 

8 

Polyakov 

[29] 

19 

19 

Ghajar  &  Asadi 

[30] 

1 

8 

Dittus-Boelter  ( n  =  1/3) 

[14] 

-18 

18 

see  above 

-7 

9 

Gnielinski’s  model  in  combination  with  Haaland’ s  friction  factor  as  well  as  the  model  proposed 
by  Ghajar  &  Asadi  predict  very  well  the  super-critical  heat  transfer  coefficient  of  CO2  under 
cooling  mode  in  the  temperature  region  investigated.  A  mean  deviation  of  only  4%  was 
calculated  for  Gnielinski’s  model. 

The  fit  of  the  rather  complicated  model  proposed  by  Polyakov  was  worse  than  the  values 
according  to  the  standard  Dittus-Boelter  correlation.  However,  the  fit  of  the  Dittus-Boelter 
correlation  using  a  Prandtl  exponent  of  n  =  0.4  (<Jmea„  =  9%)  was  much  better  than  when  using 
n  =  1/3  (amean  =  18%). 


2  The  Ghajar  and  Asadi  correlation  comes  with  its  own  equation  of  state.  This  has  not  been  implemented  here.  The 
IUPAC  equations  [42,43]  have  been  used  as  for  all  calculations.  It  might  be  that  the  Ghajar  and  Asadi  correlation 
would  have  fit  better  to  the  test  data  if  its  “own”  equation  of  state  had  been  used. 
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Comparison  with  Gnielinski  (VDI) 


/ 
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Comparison  with  Gnielinski 
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Comparison  with  simple  Gnielinski 
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+  p90  m600  qlO  xp90  m1200  q20  Dp90  m900  q20  Op80  m900  q20  Op9100  m900  q20  Ap90  m600  q20 
Figure  18:  Comparison  of  measured  heat  transfer  coefficients  with  calculation  models 
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As  one  can  see  from  Figure  18,  most  of  the  calculated  values  fitted  the  experimental  data  within 
a  deviation  of  ±15%.  Only  the  Polyakov  correlation  overestimated  the  heat  transfer  coefficient 
(Figure  19)  while  the  Dittus-Boelter  correlation  underestimated  it  (Figure  20). 

The  following  subsections  briefly  describe  the  qualitative  dependence  of  the  heat  transfer 
coefficient  on  the  temperature  and  analyse  the  influence  of  heat  flux,  mass  flow  rate,  and 
pressure.  Figure  21  to  Figure  23  show  a  comparison  of  the  measured  values  with  Gnielinski’s 
correlation  (curves).  In  the  figures,  the  mean  temperatures  of  each  experiment  is  indicated  with 
markers  while  the  bars  indicate  the  “overall  cooling  processes”  (inlet  temperature  to  outlet 
temperature). 


Comparison  with  Polyakov 
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Figure  19:  Comparison  of  measured  heat  transfer  coefficients  with  calculation  models 
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Comparison  with  Dittus-Boelter(n  =  1/3) 


+15%  0% 
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Comparison  with  Dittus-Boelter  (n  =  0.4) 
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Figure  20:  Comparison  of  measured  heat  transfer  coefficients  with  calculation  models 

Influence  of  heat  flux: 

The  influence  of  the  heat  flux  on  the  heat  transfer  coefficient  was  rather  small,  as  can  be  seen 
from  Figure  21.  Gnielinski’s  correlation  for  the  Nusselt  number  takes  into  account  the  heat  flux 
by  means  of  the  factor  (Pr/Prw)011.  Since  the  wall  temperature  changes  with  heat  flux,  the 
Prandtl  number  at  wall  temperature  Prw  changes  as  well.  However,  due  to  the  small  exponent  of 
0.1 1,  the  correction  factor  does  not  vary  much. 

If  the  heat  flux  rises,  the  heat  transfer  coefficient  makes  an  almost  parallel  translation  to  lower 
temperatures.  This  means  that  below  the  pseudo-critical  temperature  -  indicated  with  the  peak  of 
the  heat  transfer  coefficient  -  the  heat  transfer  coefficient  rises  slightly,  and  above  it  becomes 
slightly  lower  in  case  of  a  higher  heat  flux.  A  closer  investigation  showed  that  if  the  heat  flux  is 
increased  from  10,000  to  20,000  W/m2  the  (calculated)  variation  of  the  heat  transfer  coefficient 
corresponds  to  an  average  change  of  1.4%  in  the  temperature  region  from  10  to  80°C. 
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a  p91  m600  qlO  ♦  p91  m600  q20 


Figure  21:  Influence  of  the  heat  flux  on  the  heat  transfer  coefficient 
The  comparison  of  measurement  and  calculation  showed  a  high  accuracy  for  both  series. 


Influence  of  mass  flow  rate: 

In  contrast  to  that  of  the  heat  flux,  the  influence  of  the  mass  flow  rate  was  significant.  A  higher 
mass  flow  rate  of  course  lead  to  a  higher  heat  transfer  coefficient.  This  can  be  seen  from  Figure 
22.  The  calculations  showed  that  the  average  heat  transfer  coefficient  rose  by  about  42%  if  the 
mass  flux  was  increased  from  600  to  900  kg/(m2-s),  and  by  82%  if  the  mass  flow  was  increased 
to  1200  kg/(m2-s). 

The  agreement  between  calculation  and  measurement  was  quite  satisfactory  (Figure  22).  Only  at 
high  temperatures  could  a  moderate  deviation  be  seen. 


A  p91  m600  q20  ♦  p91  m900  q20  ■  p91  ml  200  q20 


Figure  22:  Influence  of  the  mass  flow  on  the  heat  transfer  coefficient 
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Influence  of  pressure: 

The  change  of  the  pseudo-critical  temperature  with  pressure  is  of  special  interest  (Table  3).  Since 
the  heat  transfer  coefficient  reaches  a  maximum  at  the  pseudo-critical  temperature,  this 
maximum  will  change  with  varying  pressure. 

Figure  23  shows  this  fact.  The  maximum  heat  transfer  increases  the  closer  the  pressure  comes  to 
the  critical  pressure  of  the  refrigerant  ( Pcrmcai  =  73.83  bar).  At  81  bar  the  maximum  value  is 
higher  than  21,000  W/(m2-K)  and  at  101  bar  the  maximum  is  about  10,000  W/(m2-K)  at  the 
pseudo-critical  temperature.  This  corresponds  to  a  rise  of  about  68%  if  the  pressure  decreases 
from  91  to  81  bar.  If  the  pressure  rises  from  91  to  101  bar,  the  maximum  heat  transfer  coefficient 
decreases  by  about  21%.  At  very  low  and  very  high  temperatures,  i.e.  away  from  the  pseudo- 
critical  temperature,  the  influence  of  pressure  is  small.  This  can  be  seen  by  the  asymptotic 
behaviour  of  the  curves  below  about  25°C  and  above  80°C. 
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Figure  23:  Influence  of  the  pressure  on  the  heat  transfer  coefficient 


Again,  the  deviation  between  measurement  and  calculation  was  small. 

5.1.3  Pressure  drop  at  super-critical  pressures 

Similarly  to  what  was  done  in  the  previous  subsection,  here  the  measured  pressure  drop  data 
were  compared  with  calculation  models  for  the  frictional  pressure  drop3.  The  selected  models 
are: 


•  Blasius  [28]:  This  is  the  common  correlation  for  estimating  the  pressure  drop  in  a  smooth 
tube  and  turbulent  flow. 

•  Colebrook  &  White  [28]:  This  implicit  form  for  the  friction  factor  gives  values  in  accordance 
to  the  well  known  Moody-Diagram.  It  can  be  applied  in  the  transition  region  between  non- 
turbulent  and  turbulent  flow. 


3  The  pressure  change  due  to  deceleration  is  neglected  here,  because  calculations  have  shown  that  this  is  low 
compared  to  the  frictional  pressure  drop. 
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•  Swamee  et  al.  have  published  an  explicit  form  of  the  Colebrook  &  White  friction  factor  [31]. 

(Remark:  For  the  latter  two,  the  tube  roughness  was  assumed  to  be  1-10'6  m.) 

The  deviations  of  the  calculated  values  from  the  measured  data  are  shown  in  Table  5.  As  can  be 
seen,  the  models  predicted  the  pressure  drop  quite  accurately.  The  mean  deviation  was 
approximately  2%  for  both  the  Colebrook  &  White  and  the  Swamee  model.4 


Table  5:  Deviation  of  the  measured  pressure  drop  to  calculation  models 


MB 

&mean  (%) 

,  •-  ■-  -  -  -  -  . . *  ""ws*  ^ , 

'  ’■  -  ■*,  Rif, cine:  ! 

v  W  .  *  f  * 

[28] 

-6 

6 

: v- wzfsf:  ■<■... 

Coolebrook  &  White 

[28] 

1 

2 

Swamee 

[31] 

1 

2 

Only  the  Blasius  correlation  slightly  underestimated  the  pressure  drop.  But  this  could  be 
expected,  because  the  present  MPE  tube  made  of  aluminium  had  a  certain  roughness. 

Figure  24  visualises  the  good  fit  of  the  models  to  the  experiments.  Almost  every  measured  point 
is  located  on  the  line  representing  a  deviation  of  0%.  Figure  25  to  Figure  27  show  the  influence 
of  heat  flux,  mass  flow  rate,  and  absolute  pressure  on  the  pressure  drop. 

As  can  be  seen  from  Figure  25,  there  was  no  noticeable  influence  on  the  pressure  drop  of  the 
heat  flux.  In  contrary  to  that,  the  mass  flow  rate  had  a  considerable  effect  (Figure  26).  As  stated 
above,  the  heat  transfer  coefficient  rose  by  about  42%  if  the  mass  flow  was  changed  from  600  to 
900  kg/(m2-s)  and  by  about  82%  if  it  was  increased  to  1200  kg/(m2-s).  However,  at  the  same  time 
the  pressure  drop  rose  about  100%  and  240%,  respectively. 

In  general,  the  dependence  of  the  pressure  drop  on  the  temperature  showed  a  change  near  the 
pseudo-critical  temperature.  The  steep  gradient  near  this  temperature  is  caused  by  the  rapid 
change  of  the  density.  Since  the  pressure  has  a  large  influence  on  the  pseudo-critical 
temperature,  the  characteristic  of  the  pressure  drop  changes  with  the  absolute  value  of  the 
pressure  (Figure  27).  In  the  temperature  region  from  10  to  80°C,  the  mean  pressure  drop 
increased  about  17%  if  the  pressure  was  changed  from  91  to  81  bar  and  decreased  13%  if  the 
pressure  rose  from  91  to  101  bar.  For  the  heat  transfer  coefficient  the  corresponding  changes 
were  68  and  21%,  respectively. 


4  It  should  be  noted  that  the  Swamee  and  the  Colebrook  and  White  correlations,  as  well  as  the  Haaland  correlation 
mentioned  in  the  previous  subsection,  are  expected  to  give  essentially  the  same  results,  as  they  are  all  equations 
made  to  fit  the  Moody  diagram. 
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A  p91  m600  qio  ♦  p91  m600  q20 


Figure  25:  Influence  of  the  heat  flux  on  the  pressure  drop  (A1  =  503  mm) 


A  p91  m600  q20  ♦  p91  m900  q20  ■  p91  ml 200  q20 


Figure  26:  Influence  of  the  mass  flow  on  the  pressure  drop  (A1  =  503  mm) 


Figure  27:  Influence  of  the  absolute  pressure  on  the  pressure  drop  (A1  =  503  mm) 


5.2  Subcritical  Conditions  (Evaporation) 


5.2.1  General 

Even  though  the  critical  temperature  of  C02  (31.1°C,  73.8  bar)  is  quite  low,  heat  absorption 
takes  place  at  subcritical  conditions  (evaporation)  in  heat  pump  and  refrigeration  applications. 
The  reduced  temperature  {T/Tmd  and  pressure  {p/p end  are  quite  high,  however. 

At  saturation  temperatures  between  0  and  20°C,  C02  has  a  high  vapour  density;  thus  the  density 
ratio  between  the  liquid  and  vapour  phase  is  significantly  lower  for  C02  than  for  commonly  used 
refrigerants.  Furthermore,  the  surface  tension  of  C02  in  this  important  temperature  region  is  very 
low. 

Taking  these  facts  into  account,  it  was  decided  to  analyse  the  heat  transfer  and  pressure  drop 
characteristics  of  C02  at  temperatures  between  0  and  20°C.  Other  reasons  for  concentrating  on 
this  temperature  were  focusing  on  air  conditioning  applications  and  using  water  as  coolant. 
Figure  28  shows  the  evaporation  process  in  a  temperature-entropy  diagram  (see  Figure  4  for  the 
flow  circuit  of  the  test  rig). 


Slightly  sub-cooled  liquid  entered  the  preheat  section  (1),  where  it  was  partly  evaporated  to  get 
the  desired  vapour  fraction  at  test  section  inlet  (2).  In  the  test  section,  refrigerant  partly 
evaporated  due  to  the  heat  supplied  by  the  heat  source  (water).  After  leaving  the  test  section  (3), 
refrigerant  entered  the  condenser,  where  it  was  condensed  to  liquid. 

In  order  to  determine  the  influence  of  mass  flux  (m),  the  evaporation  temperature  (tevap),  and 
heat  flux  (q )  on  the  ‘local’  heat  transfer  coefficient  (a),  tests  were  carried  out  at  different  vapour 
fractions  (xmean).  Table  6  gives  an  overview  of  the  test  conditions  and  the  experimental 
tolerances. 
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Table  6:  Test  conditions  and  experimental  tolerances  during  the  evaporation  tests 


Parameter 

min 

max 

max.  deviation 

m  (kg/(m2-s)) 

200 

600 

+6%  /  -7% 

tevap  (  C) 

0 

20 

±0.2  K 

q  (W/m2) 

5,000 

20,000*} 

+5%  /  -10% 

some  tests  at  higher  heat  fluxes. 


In  total,  88  points  were  collected.  Besides  the  systematic  analysis  of  these  data,  comparisons 
with  ‘common’  calculation  models  were  carried  out.  The  results  are  presented  in  the  following 
sections. 


5.2.2  Heat  T ransfer  at  evaporation 

In  contrast  to  the  supercritical  heat  transfer  coefficient,  no  specific  calculation  model  for 
evaporating  CO2  was  available  in  the  literature.  In  [13]  it  was  found  that  the  model  proposed  by 
VDI  [28]  fitted  well  to  the  measurements  published  in  [1].  But  these  measurements  were  carried 
out  in  a  tube  with  an  inner  diameter  of  7  mm.  For  the  present  tube  this  model  cannot  be  applied, 
because  of  the  small  diameter  of  0.787  mm. 

Six  models  were  chosen  for  the  comparison  with  the  experimental  data.  Figure  29  shows  the 
quantitative  comparison  between  all  measured  data  and  these  calculation  models,  and  Table  7 
lists  the  corresponding  deviations. 


Table  7:  Deviation  between  measured  and  calculated  heat  transfer  coefficients,  all  data 


Model 

".Reference,  | 

Slipcevic  (1987) 

[32] 

60 

68 

Kandlikar  (1990) 

[33] 

49 

55 

Kattan  (1996) 

[34] 

91 

92 

Shah  (1982) 

[35] 

4 

61 

Gungor  &  Winterton  (1986) 

[36] 

125 

56 

Gungor  &  Winterton  (1987) 

[37] 

14 

117 

As  one  can  see,  there  was  a  considerable  discrepancy  between  calculation  and  experiment. 
Especially,  in  case  of  low  measured  heat  transfer  coefficients,  the  models  delivered  too  high 
values.  If  one  analyses  the  data  in  depth,  the  major  discrepancy  is  seen  to  occur  at  high  mass 
fluxes.  Therefore,  Table  8  shows  the  deviations  for  the  points  measured  at  mass  fluxes  less  than 
or  equal  to  300kg/(m*s).  Then  the  remaining  40  points  show  a  better  correspondence. 
Especially  Slipcevic  and  Kandlikar  fitted  the  measurements  very  well  at  low  mass  fluxes. 
However,  the  large  deviation  at  high  mass  fluxes  has  to  be  considered. 
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Figure  29:  Comparison  of  measured  heat  transfer  coefficients  with  calculation  models 
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Table  8:  Deviation  between  measured  and  calculated  heat  transfer  coefficients,  m  <  300  kg/(m2-s) 


Model 

Oavg  (%) 

Slipcevic  (1987) 

[32] 

-6 

15 

Kandlikar  (1990) 

[33] 

0 

15 

Kattan  (1996) 

[34] 

8 

22 

Shah  (1982) 

[35] 

-41 

42 

Gungor  &  Winterton  (1986) 

[36] 

44 

44 

Gungor  &  Winterton  (1987) 

[37] 

-33 

33 

Figure  30  shows  the  results  obtained  from  the  initial  evaporation  tests.  The  conditions  applied 
were:  m  =  600  kg/(m2*s),  q  =  20,000  W/m2,  and  different  evaporation  temperatures.  As  can  be 
seen,  the  heat  transfer  coefficient  dropped  significantly  at  high  vapour  fractions  (x).  At  low 
vapour  fractions  it  was  about  13,000  W/(m2-K)  at  0°C,  15,000  W/(m2-K)  at  10°C,  and 
22,000  W/(m2-K)  at  20°C  evaporation  temperature.  At  vapour  fractions  higher  than  0.55,  the 
measured  values  were  lower  than  5,500  W/(m2-K)  for  all  evaporation  temperatures. 


_  -O — t  0/q20/m600  — * — H0/q20/m600  ---a---  t20/q20/m600 


Figure  30:  Measured  heat  transfer  coefficients  depending  on  the  vapour  fraction  at  a  high  mass 

flux  (one  refrigerant  inlet) 

Initially,  the  test  section  was  supplied  with  CO2  from  one  side  of  the  manifold,  i.e.  an 
asymmetric  flow.  This  was  believed  to  cause  a  separation  of  the  liquid  and  vapour  phase  in  the 
manifold,  which  could  cause  non-uniform  supply  of  liquid  to  the  parallel  ports  of  the  test  tube. 
This  is  shown  schematically  in  Figure  31,  which  shows  a  possible  solution  as  well. 

It  was  decided  to  redesign  the  manifolds  so  that  the  fluid  could  enter  and  leave  the  test  section 
symmetrically,  i.e.  through  two  inlets  and  two  outlets.  After  the  reconstruction  of  the  manifolds, 
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the  previous  test  was  repeated  at  an  evaporation  temperature  of  10°C.  Figure  32  shows  the  result 
for  the  ‘new’  test  section  (indicated  with  ‘2  inlets’)  in  comparison  with  the  ‘old’  results. 


C02  -  vapour/liquid  C02  -  vapour/liquid 


C02  -  vapour/liquid 

a)  b) 

Figure  31:  Possible  problem  and  suggested  solution 


* —  tl  0/q20/m600  ---a---  tl  0/q20/m600/2inlets 


Figure  32:  Heat  transfer  coefficient  depending  on  the  vapour  fraction  at  a  high  mass  flux  (one  vs. 

two  refrigerant  inlets) 

As  can  be  seen  from  this  figure,  there  was  more  or  less  no  change.  Thus  the  conclusion  was  that 
the  asymmetric  design  of  the  inlet  gave  uniform  liquid  supply  to  the  ports  and  the  reason  for  the 
low  heat  transfer  coefficient  at  high  vapour  fractions  was  yet  to  be  found. 

It  was  decided  to  run  the  same  test  once  more,  but  this  time  at  a  reduced  heat  flux  of 
5,000  W/m2.  The  result  did  not  change  significantly  at  high  vapour  fractions,  however.  To  find 
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out  if  the  high  mass  flux  of  600  kg/(m2-s)  caused  this  characteristic,  a  new  test  series  was  carried 
out. 

During  these  tests,  the  vapour  fractions  at  the  inlet  and  outlet  of  the  test  section  were  kept 
constant,  while  the  mass  flux  was  changed.  This  was  realised  by  varying  the  heat  flux.  The  heat 
fluxes  required  to  keep  the  vapour  fractions  constant  were  10,000  to  26,000  W/m2  for  the  series 
with  an  inlet  vapour  fraction  of  xin  =  0.3  and  an  outlet  vapour  fraction  of  xout  =  0.92,  and  14,000 
to  42,000  W/m2  for  the  series  with  xin  =  0  and  xoa,  =  0.9. 

Figure  33  shows  the  result.  As  can  be  seen,  the  influence  of  the  mass  flux  was  significant.  At 
both  tests,  an  increase  from  200  to  300  kg/(m2-s)  resulted  in  a  higher  heat  transfer  coefficient.  At 
higher  velocities,  a  significant  decrease  was  observed. 


m  [kg/(m2  s)] 


Figure  33:  Heat  transfer  coefficient  depending  on  the  mass  flux  (constant  vapour  fraction  at  the 

inlet  and  at  the  outlet) 


The  higher  heat  transfer  coefficient  in  case  of  the  lower  inlet  vapour  fraction  can  be  explained  by 
the  effect  of  the  heat  flux  (see  below). 


A  possible  explanation  for  the  decreasing  heat  transfer  coefficient  with  increasing  mass  flux  is 
the  dry-out  phenomenon,  in  which  the  flow  itself  causes  the  liquid  which  is  wetting  the  wall  to 
be  entrained  in  the  core  of  the  flow  at  high  mass  fluxes.  VDI  [28]  recommends  a  rather 
complicated  calculation  model  for  the  so-called  critical  vapour  fraction  {xcrmcai)  where  dry-out 
occurs.  This  model  is  mainly  based  on  experiments  with  water.  It  is  required  to  ‘translate’  the 
CO2  properties  into  equivalent  water  properties  so  that  the  calculation  can  be  carried  out.  This 
was  the  main  reason  why  it  was  unclear  if  this  model  could  be  applied  for  CO2. 

Figure  34  shows  the  result  of  the  calculations  for  different  evaporation  temperatures  (0,  10,  and 
20°C)  and  heat  fluxes  (10,000  and  20,000  W/m2).  In  case  of  20°C  saturation  temperature  and 
moderate  mass  fluxes,  the  figure  shows  three  lines  for  the  critical  vapour  fraction.  In  this  mass 
flux  range,  gravity  is  expected  to  influence  the  dry-out  phenomenon.  The  upper  dashed  line  in 


i 

j 
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the  figure  represents  the  bottom  and  the  lower  line  the  top  of  the  tube.  The  middle  line  gives  the 
average  value.  In  other  words,  dry-out  at  the  top  occurs  earlier  than  at  the  bottom  of  the  tube. 


Figure  34:  Critical  vapour  fraction  according  to  VDI  [28] 


This  figure  confirms  the  speculation  that  dry-out  on  the  heat  transferring  surface  can  occur  even 
at  rather  low  vapour  fractions.  For  a  mass  flux  of  600  kg/(m2-s)  and  a  heat  flux  of  20,000  W/m2, 
this  model  predicts  dry-out  at  a  vapour  fraction  of  0.77  at  0°C,  0.6  at  10°C,  and  0.47  at  20°C 
evaporation  temperature.  At  a  mass  flux  of  300  kg/(m2-s)  the  corresponding  values  are  0.96, 
0.76,  and  0.6,  respectively.  As  may  be  observed  by  comparing  the  first  diagram  with  the  second 
diagram  of  Figure  34,  the  lower  the  heat  flux  is,  the  higher  the  critical  vapour  fraction  becomes. 

The  measurements  shown  in  Figure  32  completely  confirmed  the  calculation.  As  the  vapour 
fraction  became  higher  than  0.6,  a  significant  decrease  of  the  heat  transfer  coefficient  occurred. 
At  an  evaporation  temperature  of  0°C  the  model  probably  overestimates  the  critical  vapour 
fraction,  as  can  be  seen  from  a  comparison  with  Figure  30. 

Since  dry-out  is  a  quite  important  phenomenon  with  respect  to  heat  exchanger  design,  more 
measurements  were  needed. 
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Influence  of  mass  flux: 

Tests  were  carried  out  at  an  evaporation  temperature  of  0°C,  10°C,  and  20°C  and  at  mass  fluxes 
of  200,  300,  and  400  kg/(m2-s).  At  all  test  runs  the  heat  flux  was  kept  constant  at  10,000  W/m2 
Figure  35  shows  the  results. 

From  the  first  diagram  (0°C)  one  can  see  that  at  low  vapour  fractions,  neither  the  mass  flux  nor 
the  vapour  fraction  had  any  influence  on  the  heat  transfer  coefficient.  At  all  velocities,  the 
measured  heat  transfer  coefficient  was  about  9000  W/(m2K).  But  again,  at  high  vapour  fractions 
and  a  mass  flux  of  400  kg/(m2-s),  the  heat  transfer  coefficient  was  significantly  reduced. 

A  similar  behaviour  can  be  seen  at  10°C  (second  diagram);  at  low  vapour  fractions  the  heat 
transfer  coefficient  was  about  10,500  W/(m2-K)  and  independent  of  the  mass  flux.  In  case  of 
200  kg/(m  -s)  this  value  was  measured  at  all  vapour  fractions.  But  already  at  300  kg/(m  -s)  and 
high  vapour  fractions,  a  slight  drop  can  be  seen.  At  20°C  (third  diagram)  and  low  vapour 
fractions  a  value  of  about  15,000  W/(m2-K)  was  measured.  The  drop,  however,  occurred  at  lower 
vapour  fractions.  These  results  are  -  except  for  the  measurement  at  20°C  /  200  kg/(m2-s)  - 
completely  in  accordance  with  the  calculated  results  for  the  critical  vapour  fraction  presented  in 
Figure  34. 

Influence  of  heat  flux: 

The  influence  of  heat  flux  on  the  heat  transfer  characteristic  of  C02  is  of  high  interest,  as  already 
described  in  [1],  Within  this  project,  the  influence  of  heat  fluxes  ranging  from  5000  to 
15,000  W/m2  was  investigated.  During  these  experiments,  the  mass  flux  was  kept  constant  at 
300  kg/(m2-s)  while  the  saturation  temperature  was  varied.  In  Figure  36  the  experimental  results 
are  shown. 

At  0°C  saturation  temperature  (first  diagram),  a  slightly  higher  heat  transfer  coefficient  was 
measured  when  the  heat  flux  was  high.  At  10°C  (second  diagram),  the  enhancement  due  to  a 
higher  heat  flux  was  larger.  At  high  vapour  fractions  an  exact  analysis  was  not  possible  because 
of  the  different  range  of  vapour  fraction  of  each  experiment  depending  on  the  heat  flux.  At  20°C, 
the  change  in  vapour  fraction  is  large  in  case  of  a  heat  flux  of  15,000  W/m2.  However,  the 
decreased  heat  transfer  coefficient  at  a  high  vapour  fraction  may  confirm  the  fact  that  due  to  the 
higher  heat  flux,  dry-out  occurs  at  a  lower  vapour  fraction.  This  would  mean  that  the 
enhancement  due  to  an  increased  heat  flux  does  not  compensate  for  the  drop  due  to  dry-out 
occurring  at  a  lower  vapour  fraction. 
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Figure  35:  Influence  of  the  mass  flux  on  the  heat  transfer  coefficient  at  different  evaporation 

temperatures 


(m2  K) 


i —  t20/q1 0/m300  — a-  ■  -  t20/q15/m300 

Figure  36:  Influence  of  the  heat  flux  on  the  heat  transfer  coefficient  at  different  evaporation 

temperatures 


Influence  of  evaporation  temperature: 

The  influence  of  evaporation  temperature  on  the  heat  transfer  characteristics  of  C02  was  already 
mentioned.  Figure  37  compares  measured  heat  transfer  coefficients  at  constant  heat  flux  and 
mass  flux,  and  varying  saturation  temperature.  As  may  be  observed,  the  heat  transfer  coefficient 
was  significantly  higher  for  the  highest  saturation  temperature  at  low  vapour  fractions.  While  the 
heat  transfer  coefficient  was  about  9,000  W/(m2-K)  at  0°C,  it  was  higher  than  15,000  W/(m2  K) 
at  20°C.  At  high  vapour  fractions,  the  situation  is  reversed  because  of  the  significant  drop  at  high 
temperatures. 


Figure  37:  Influence  of  evaporation  temperature  on  the  heat  transfer  coefficient 

5.2.3  Pressure  Drop  at  evaporation 

For  the  comparisons  of  measured  data  to  calculated  values,  the  following  correlations  were  used: 

•  Thome  (1997):  Thome  [38]  recommended  the  Friedel  correlation  using  the  Premoli  approach 

for  the  void  fraction  to  calculate  the  frictional  and  acceleration  pressure  drop  at 
evaporation.  As  shown  in  [13],  this  model  fits  well  to  measured  data  carried  out  in 
a  7  mm  round  tube  (published  by  Bredesen  at  al  in  [1]). 

•  Fuchs&Neraas  (1993):  Neraas  published  this  model  in  [39].  It  is  based  on  the  work  of  Fuchs 

[40].  Bredesen  at  al.  [1]  proposed  this  model  for  the  estimation  of  the  frictional 
pressure  drop  of  evaporating  C02. 

•  VDI  (1994):  In  [28]  a  rather  sophisticated  model  for  the  prediction  of  the  pressure  drop  can 

be  found. 

•  Single-phase  model:  This  model  uses  a  homogeneous  approach  by  means  of  the  mean 

density  for  the  two-phase  flow.  The  frictional  pressure  drop  is  then  calculated  for 
the  fieri ve  single-phase  flow  according  to  the  Colebrook&White  correlation  [28], 
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Table  9  shows  the  resulting  deviations,  and  Figure  38  visualises  the  fit  of  the  models  to  the 
measurements.  Thome’s  recommendation  gave  the  best  results.  The  mean  deviation  is 
nevertheless  22%.  All  other  models  -  and  especially  the  Fuchs&Neraas  model  that  was  found  to 
give  good  correlation  to  measured  data  for  the  7  mm  tube  -  underestimated  the  pressure  drop 
significantly. 


Table  9:  Deviation  between  measured  and  calculated  pressure  drops 


Model 

Reference 

■tg 

i 

1 

Thome  (1997) 

[38] 

-7 

22 

Fuchs&Neraas  (1993) 

[39] 

-56 

57 

VDI  (1994) 

[28] 

-41 

43 

Single-phase  model 

see  above 
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Figure  38:  Comparison  of  measured  pressure  drop  data  with  calculation  models 

In  the  following  figures,  some  of  the  measured  data  are  compared  with  the  results  calculated 
according  to  Thome  [38]. 
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The  influence  of  mass  flux  on  the  pressure  drop  at  a  saturation  temperature  of  10°C  can  be  seen 
in  Figure  39.  While  the  measurements  at  200  kg/(m2-s)  showed  a  moderate  pressure  drop  at  all 
vapour  fractions,  the  pressure  drop  at  400kg/(m2-s)  increased  significantly  with  the  vapour 
fraction.  As  one  can  see,  the  calculated  values  did  not  correlate  very  well  to  the  experimental 
data.  Especially  at  low  vapour  fractions  the  deviation  was  high. 

Finally,  Figure  40  shows  the  influence  of  the  saturation  temperature  on  the  pressure  drop.  The 
mass  flux  was  constant  at  300 kg/(m2s)  during  these  tests.  With  decreasing  saturation 
temperature,  the  vapour  volume  increases,  and  this  is  the  main  reason  for  the  remarkably  higher 
pressure  drop  at  0°C  compared  to  that  at  20°C.  Regarding  the  deviation  between  measurement 
and  calculation,  one  can  see  that  the  model  underestimates  the  pressure  drop  in  case  of  low 
temperatures. 
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Figure  39:  Influence  of  the  mass  flux  on  the  two-phase  pressure  drop  (A1  =  503  mm) 
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Figure  40:  Influence  of  the  saturation  temperature  on  the  two-phase  pressure  drop  (A1  =  503  mm) 
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6  ASSESSMENT  OF  RESULTS 


In  this  chapter,  the  evaporation  heat  transfer  data  presented  in  the  previous  chapter  will  briefly  be 
assessed  regarding  their  consequences  for  the  design  of  C02-based  air  conditioning  systems.  It 
will  be  done  by  means  of  comparing  previous  capacity  measurements  of  evaporators  conducted 
in  the  laboratory  at  SINTEF  Energy  Research  with  calculations  of  evaporators  using  correlations 
and  the  heat  transfer  data  obtained  in  this  work. 

6.1  Previous  evaporator  measurements 

In  a  previous  project,  the  performance  of  vehicle  air  conditioning  evaporators  were  measured. 
The  prototype  evaporator  consisted  of  extruded  aluminium  multiport  tubes  of  similar  geometry 
to  the  one  used  in  the  present  work.  A  photograph  of  the  evaporator  is  shown  in  Figure  41. 


Figure  41:  Prototype  automotive  aluminium  multiport  extruded  tube  evaporator. 

The  evaporator  inlet  air  temperature  ranged  from  26-48  °C,  the  inlet  relative  air  humidity  was 
varied  between  29  and  28  %  and  the  air  flowrate  was  between  386  and  394  Sm3/h. 

The  C02  mass  flux  ranged  from  270  to  900  kg/(m2s)  (it  also  varied  internally  in  the  evaporator 
due  to  the  varying  number  of  parallel  tubes).  The  evaporation  temperature  varied  from  —1  to  15 
°C  and  the  heat  flux  ranged  from  3,500  and  10,000  W/m2. 
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6.2  Evaporator  simulations 

The  computer  program  HXSIM  was  developed  at  SINTEF  Energy  Research  to  assist  the  design 
of  compact  evaporators,  gas  coolers  and  condensers.  It  simulates  the  heat  exchangers  considering 
a  great  range  of  parameters. 

Simulations  of  the  test  cases  mentioned  in  the  previous  section  were  performed  with  three 
different  calculation  methods  for  the  CCVside  heat  transfer  coefficient: 

1.  The  C02-side  heat  transfer  coefficient  was  calculated  with  the  simplified  Gungor  and 
Winterton  (1987)  correlation  [37] .  (The  Kandlikar  (1990)  correlation  was  not  implemented  in 
HXSIM). 

2.  The  CC>2-side  heat  transfer  coefficient  was  calculated  with  the  simplified  Gungor  and 
Winterton  (1987)  correlation,  but  the  drop  in  heat  transfer  coefficient  due  to  dry-out  was 
accounted  for  by  the  following  modification:  For  vapour  fractions  larger  than  the  critical 
vapour  fraction  given  by  VDI  [28],  the  nucleation  boiling  factor  (boiling  number)  was  set  to 
zero.  Then  the  effective  heat  transfer  coefficient  was  calculated  as  outlined  in  [41]. 

3.  The  CO2  -side  heat  transfer  coefficient  was  based  on  the  experimental  data  presented  in  this 
work.  For  each  of  the  evaporator  test  points,  a  heat  transfer  curve  for  a  suitable  evaporation 
temperature  and  mass  flux  was  found.  Then  the  heat  transfer  coefficient  was  given  as  a 
function  of  vapour  fraction  to  the  HXSIM  program.  The  HXSIM  program  thus  used  one¬ 
dimensional  interpolation  of  the  test  data  to  find  the  CC>2-side  heat  transfer  coefficient  when 
simulating  the  evaporator. 


Figure  42:  Measured  and  calculated  cooling  capacities  of  MPE  evaporator. 
The  symbols  are  explained  in  the  text 
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The  results  of  the  calculations  are  shown  in  Figure  42.  The  capacities  calculated  using  the 
Gungor  and  Winterton  (1987)  correlation  (labelled  GW87)  all  lay  above  the  measured  values. 
The  thick  line  is  a  linear  curve-fit  to  these  values.  At  6,000  W  capacity,  GW87  yielded  a  capacity 
17  %  higher  than  the  measured  one.  GW87_cr  denotes  the  Gungor  and  Winterton  (1987) 
correlation  modified  to  take  account  for  the  drop  in  the  heat  transfer  coefficient  caused  by  dry¬ 
out.  These  values  lay  only  slightly  lower  than  the  GW87  values.  The  capacities  obtained  using 
the  measured  heat  transfer  data  are  labelled  MPE.  The  dotted  line  is  a  linear  curve-fit  through 
these  points.  At  6,000  W,  the  MPE  data  yielded  a  capacity  14  %  higher  than  measured,  i.e.  3 
percentage  points  lower  than  what  was  given  by  GW87. 

6.3  Discussion 

The  evaporator  cooling  capacities  obtained  using  measured  MPE  heat  transfer  data  were  only 
slightly  lower  than  those  obtained  using  the  Gungor  and  Winterton  (1987)  correlation.  This  is 
perhaps  surprising,  having  seen  in  the  previous  chapter  that  the  average  deviation  between  this 
correlation  and  the  measured  heat  transfer  coefficients  was  14  %,  and  that  the  mean  deviation 
was  as  high  as  1 17  %.  Generally  one  may  assume  the  possible  discrepancy  between  calculated 
and  measured  heat  transfer  coefficients  to  be  in  the  range  of  20  %. 

The  most  important  reason  for  the  small  differences  between  the  calculated  evaporator  capacities 
for  different  C02-side  heat  transfer  coefficients,  was  that  as  much  as  around  80  %  of  the  heat 
transfer  resistance  was  on  the  air-side.  Thus  one  can  roughly  say  that  a  20  %  increase  in  the  C02- 
side  heat  transfer  coefficient  only  results  in  a  4  %  increase  in  the  transferred  heat. 

The  most  important  cause  for  most  of  the  heat  transfer  resistance  being  on  the  air-side  was  that 
the  heat  transfer  coefficient  on  the  C02-side  was  very  high,  together  with  the  fact  that  the  ratio 
between  the  air-side  surface  and  the  C02-side  surface  was  quite  low.  A  different  refrigerant  and / 
or  different  operating  conditions  leading  to  a  lower  refrigerant-side  heat  transfer  coefficient 
would,  therefore,  result  in  variations  in  the  refrigerant-side  heat  transfer  coefficient  having  a 
greater  influence  on  the  evaporator  cooling  capacity  than  seen  in  this  comparison. 

Important  reasons  for  the  about  15  %  deviation  between  simulated  and  measured  evaporator 
performance  include: 

•  Effects  of  retained  water  on  the  air-side. 

•  Presence  of  lubricant  in  the  refrigerant  flow. 

The  evaporator  discussed  here  was  a  prototype  with  a  rather  small  air-side  area  compared  to  that 
on  the  C02-side  (Aai/ACo2  =  6).  Later  evaporators  would  probably  be  constructed  with  a  larger 
air-side  area,  thus  approaching  an  equal  distribution  of  the  heat  transfer  resistances  between  the 
two  sides.  In  such  an  instance,  a  20  %  increase  in  the  C02-side  heat  transfer  coefficient  would 
lead  to  a  roughly  10  %  increase  in  the  transferred  heat,  which  is  not  insignificant.  The  same 
would  be  the  case  for  e.g.  a  water-heated  evaporator,  where  heat  transfer  coefficients  on  both 
sides  would  be  high. 
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7  CONCLUSIONS 


Within  this  project,  the  heat  transfer  and  pressure  drop  characteristic  of  C02  (R-744)  was 
investigated.  The  focus  was  on  the  heat  rejection  process  at  supercritical  pressures  and  the  heat 
absorption  process  at  subcritical  pressures  (evaporation).  The  test  tube  was  a  so-called  MultiPort 
Extruded  tube  with  micro-channels.  It  had  25  parallel  ports  on  the  refrigerant-side  with  a  port 
diameter  of  0.787  mm. 

Because  of  the  low  critical  temperature  of  C02  (31°C),  the  so-called  transcritical  vapour 
compression  cycle  will  be  used  in  several  applications,  including  U.S.  Army  Environmental 
Control  Units.  Therefore,  the  test  conditions  were  selected  in  order  to  collect  knowledge  on  near- 
critical  processes. 

In  the  transcritical  cycle,  refrigerant  is  cooled  at  supercritical  pressure,  where  the  influence  of  the 
critical  point  on  the  properties  is  quite  large,  and  it  evaporates  at  temperatures  close  to  the  critical 
temperature,  where  the  high  saturation  pressure  results  in  high  vapour  densities.  Furthermore,  the 
surface  tension  of  C02  is  very  low.  These  factors  lead  to  conditions  in  C02  equipment  differing 
considerably  from  those  of  the  hitherto  commonly  used  refrigerants. 

The  experimental  results  confirm  that  C02  offers  high  heat  transfer  coefficients  at  super-critical 
pressures.  A  comparison  of  these  data  with  common  correlations  showed  good  correspondence. 
Especially  the  Nusselt  number  based  on  the  Gnielinski  correlation  in  combination  with 
Haaland’s  friction  factor  correlated  very  well  with  the  experimental  data.  This  model  can  be 
found  in  the  Appendix. 

The  pressure  drop  of  C02  to  be  expected  in  a  refrigerant  cooler  is  rather  high,  but  due  to  the  high 
pressure  level,  the  effect  on  temperature  loss  is  moderate  [13].  A  comparison  of  the  experimental 
data  with  calculations  using  the  Colebrook  &  White  correlation  showed  satisfactory  agreement. 
This  model  is  also  outlined  in  the  Appendix. 

The  heat  transfer  characteristics  of  C02  during  evaporation  is  not  as  clear  as  the  situation  at 
supercritical  pressures.  The  effect  of  the  low  density  ratio  between  liquid  and  vapour  may 
significantly  affect  the  two-phase  flow.  In  case  of  high  mass  fluxes,  the  experiments  showed  a 
strongly  decreasing  heat  transfer  coefficient  from  a  certain  vapour  fraction  upwards.  In  the 
temperature  region  investigated,  the  ‘critical’  mass  flux  was  about  300  kg/(m2-s),  with  somw 
variation  depending  on  the  saturation  temperature  and  the  heat  flux.  A  comparison  with  a 
calculation  model  for  dry-out  showed  that  this  drop  had  to  be  expected.  The  calculated  values  for 
the  critical  vapour  fraction  corresponded  well  to  the  observations.  This  very  important  fact  has  to 
be  considered  when  designing  evaporators. 
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None  of  the  investigated  heat  transfer  calculation  models  takes  this  phenomenon  into  account.  If 
the  measurements  at  high  mass  flow  rates  were  not  included,  the  models  proposed  by  Slipcevic 
and  by  Kandlikar  yielded  sufficient  accuracy.  Kandlikar’s  model  is  given  in  the  Appendix. 
However,  regarding  the  qualitative  accuracy  of  this  model,  one  has  to  remember  that  only 
measured  mean  values  have  been  compared  with  the  calculation  results. 

The  two-phase  pressure  drop  correlations  yielded  too  low  values  in  general.  The  calculation 
model  proposed  by  Thome,  the  Friedel  correlation  for  the  frictional  pressure  drop,  and  Premoli’s 
approach  for  the  void  fraction  (see  Appendix),  can  be  used  for  a  rough  estimation  of  the  pressure 
drop  to  be  expected.  But  at  low  temperatures  and  low  vapour  fractions,  the  calculated  values 
were  significantly  lower  than  the  measured  ones.  A  model  used  for  calculating  the  pressure  drop 
in  a  7  mm  round  tube  could  not  be  used  for  predicting  the  pressure  drop  in  the  micro-channel 
tube  that  was  investigated  within  this  project. 

A  variation  in  the  C02-side  heat  transfer  coefficient  had  only  a  minor  influence  on  the  evaporator 
capacity  of  a  prototype  automobile  air  conditioning  evaporator  considered  in  this  work.  On  this 
evaporator,  most  of  the  heat  transfer  resistance  lay  on  the  air-side.  For  evaporators  with  a  more 
evenly  distributed  heat  transfer  resistance  between  the  C02-side  and  ‘the  other’  side,  however,  a 
change  in  the  C02-side  heat  transfer  coefficient  will  have  an  effect  on  the  evaporator  capacity. 


All  the  results  achieved  within  this  project  are  important  for  the  design  of  C02  heat  exchangers. 
Especially  the  influence  of  the  mass  flux  and  dry-out  on  the  two-phase  heat  transfer 
characteristics  has  to  be  taken  into  consideration  in  further  work. 
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NOMENCLATURE 


Symbols 

Description 

Formula 

Unit 

A 

Area,  surface 

- 

VO 

A 

•  Wilson  Plot  constant  (slope) 

- 

- 

B 

Wilson  Plot  constant  (intersection) 

- 

- 

Bo 

Boiling  number 

- 

C 

Constant 

- 

- 

Cc 

Coefficient  of  contraction 

- 

- 

Co 

Convective  number 

- 

Cp 

Specific  isobaric  heat  capacity 

- 

J/(kgK) 

D 

Log.  Wilson  Plot  constant  (slope) 

- 

- 

dhyd 

Hydraulic  diameter 

- 

m 

E 

Log.  Wilson  Plot  constant  (intersection) 

- 

- 

F,  G,  H 

Factor 

Fa 

Property  factor 

Fa  =  0.023 -Pr 04 -/l/V8 

- 

f 

Friction  factor 

- 

- 

Fr 

Froude  number 

g 

Gravity 

- 

m/s2 

h 

Specific  enthalpy 

- 

J/kg 

h' 

Specific  enthalpy  of  saturated  liquid 

- 

J/kg 

h' 

Specific  enthalpy  of  saturated  gas 

- 

J/kg 

k 

Overall  heat  transfer  coefficient 

- 

W/(m2-K) 

Kc 

Entrance  pressure  loss  coefficient 

- 

- 

Ke 

Exit  pressure  loss  coefficient 

- 

- 

m 

Reynolds  exponent 

-  . 

- 

m 

Specific  mass  flow,  mass  flux 

m  =  M/ A 

kg/(m2-s 

M 

Mass  flow  rate 

- 

kg/s 

n 

Prandtl  exponent 

- 

- 

n 

Number  of  cells 

- 

- 

Nu 

Nusselt  number 

Nu  =  (a-dhyd)/A 

- 

P 

Pressure 

- 

Pa,  bar 

Pr 

Prandtl  number 

Pr  =  cp  -rj/A 

- 

q 

Heat  flux 

>■55. 

II 

to- 

W/m2 

Q 

Heat  transfer  rate 

- 

W 

Re 

Reynolds  number 

Re  =  ih-dhyd/T] 

- 

Rw 

Heat  transfer  resistance 

- 

(m2-K)/W 

s 

Specific  entropy 

- 

J/(kg-K) 

t 

Temperature 

- 

°C 

67 


Symbols 

Description 

Formula 

Unit 

t 

Mean  temperature 

- 

°C 

T 

Absolute  temperature 

- 

K 

t  pseudo 

Pseudo-critical  temperature 

- 

°C 

V 

Specific  volume 

1/p 

m3/kg 

We 

Weber  number 

- 

- 

X 

Vapour  fraction 

x  =  (h  -  h')/(h"  -  h') 

- 

X, 

Independent  variable  of  Wilson  Plot 

- 

- 

x2 

Independent  variable  of  log.  Wilson  Plot 

- 

- 

Yi 

Dependent  variable  of  Wilson  Plot 

- 

- 

y2 

Dependent  variable  of  log.  Wilson  Plot 

“ 

a 

Heat  transfer  coefficient 

1 

* 

II 

$ 

W/(m2-K) 

A 

Difference 

- 

- 

£ 

Void  fraction 

- 

- 

£ 

Roughness  of  the  tube  surface 

- 

m 

0 

Two-phase  multiplier 

- 

- 

rj 

Dynamic  viscosity 

- 

Pa  •  s 

A 

Friction  factor 

- 

- 

A 

Thermal  conductivity 

- 

W/(m  •  K) 

P 

Density 

- 

kg/m3 

<7 

Ratio  of  port  area  to  front  area 

®  ~  Ap.tot / ■A fir 

- 

a 

Surface  tension 

- 

N/m 

\ 

Friction  factor 

- 

- 

Subscripts  Description 

c 

cold  fluid 

c 

contraction 

calc 

calculated 

con 

post  condenser 

con 

condenser 

cool 

cooling 

cooler 

water  cooler 

e 

expansion 

el 

electrical 

eq 

equivalent 

fl 

fluid 
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Subscripts  Description 


evap 

evaporation 

f 

friction 

fr 

front 

g 

gas 

h 

heat  exchanger 

h 

hot  fluid 

heat 

heating 

h 

homogeneous 

hyd 

hydraulic 

i 

inner 

i 

inlet 

i 

cell  i 

ips 

inlet  preheat  section 

its 

inlet  test  section 

l 

liquid 

Im 

logarithmic  mean 

m 

mean 

meas 

measured 

mps 

middle  of  preheat  section 

o 

outer 

o 

outlet 

ots 

outlet  test  section 

p 

port 

ref 

refrigerant 

spec 

specific 

tot 

total 

ts 

test  section 

w 

wall 

w 

water  flow 

w 

water 

APPENDIX :  CORRELATIONS 


A.1  Super-critical  Conditions  (Cooling) 

A.1.1  Heat  Transfer 

The  mean  Nusselt  number  for  turbulent  flow  in  a  tube  (2,300  <  Re  <  1  •  106)  can  be  calculated 
according  to  Gnielinski’s  correlation  VDI  (1994): 


— -(i?e-1000)-Pr 

NUm  =■ - r = - 

1  +  12.7  -M-(Pr2'3-l) 


1  + 


(d^ 

2/3  " 

I'J 

Eq  1 


where  l  is  the  tube  length,  and /is  the  pressure  drop  factor.  Here  the  measurements  justify  the  use 
of  Haaland’s  friction  factor  (Munkejord,  1998): 


/=  -1.8-Zog 


6-9  +(e/d 


Re 


V" 


3.7 


Eq  2 


where  £is  the  tube  roughness.  Comparsions  have  shown  a  good  correspondence  of  this 
correlation  to  the  well-known  Colebrook  &  White  correlation  (Eq  5). 


The  influence  of  the  temperature  on  the  properties  can  be  taken  into  account  by 


Nu  =  Nu„ 


s  xO.ll 

f  Pr  ' 


Pr 


for  liquids 


Eq  3 


where  Pr  is  the  Prandtl  number  at  the  average  fluid  temperature  and  Prw  is  the  Prandtl  number  of 
the  fluid  at  wall  temperature  (range  of  validity  0.1  <  Pr/Prw<  10).  Since  there  is  no  noticeable 
influence  of  temperature  on  the  Prandtl  number  in  the  case  of  gaseous  fluids,  the  influence  is 
taken  into  account  by 


Nu  =  Num 

m 


r  T  V 


....for  gaseous  fluids 


Eq  4 


where  T  is  the  average  fluid  temperature  [K]  and  Tw  is  the  average  wall  temperature  [K].  For 
cooling  (T/Tw >  1)  the  exponent  n  is  zero,  for  heating  (1  >  T/Tw>  0.5)  n  depends  on  the  fluid;  for 
CO2,  n  =  0.12. 


Remark:  In  VDI  (1994)  no  comments  concerning  the  influence  of  the  wall  temperature  on  the 
properties  at  super-critical  pressures  are  made;  in  this  work  equation  Eq  3  has  been  used  at 
super-critical  pressures. 
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A.1.2  Pressure  Drop 

For  the  single-phase  pressure  drop  VDI  (1994)  recommends  the  implicit  form  for  the  friction 
factor  according  to  Colebrook  &  White: 


1  _  .  2.51  dd 

—  =  -2- log  - =  +  -!— 

f  [Re-Jf  3.71 


Eq  5 


where  8  is  the  roughness  of  the  tube  surface. 


A.2  Sub-critical  Conditions  (Evaporation) 

A.2.1  Heat  Transfer 

Kandlikar  (1990)  has  used  for  his  equation  5246  data  points  from  24  authors  for  water,  7  CFCs, 
nitrogen  and  neon.  The  equation  is  valid  for  both  horizontal  and  vertical  tubes. 

The  heat  transfer  coefficient  is  calculated  for  convective  boiling  and  nucleate  boiling,  the  larger  one 
is  taken  as  heat  transfer  coefficient  (the  ‘greater  of  the  two’  method). 

The  heat  transfer  coefficient  is  defined  by 

a  =  afcrCoc-(25-Fr'Y‘+CrBoc'F,\  Eq6 


where  Fji  is  fluid  dependent,  C1.5  are  constants,  Fr '  the  liquid  Froude  number: 


Fr'=- 


m 


P  '-g-d 


Eq  7 


The  Boiling  number  is  defined  as 


Bo  = 


m-{h"—h') 
and  the  convection  number  as 


Eq  8 


Co  = 


1 — x 


\  0.8  f 


Eq  9 


which  in  fact  is  a  modified  Martinelli  parameter. 

The  single-phase  liquid  heat  transfer  coefficient  GCi  is  calculated  using  the  equation  of  Dittus-Boelter 
a,  =  0.023  •  Ref- 8  •  Pr'°  4  ■  X7  d  Eq  10 


Interesting  is  the  fluid-dependent  parameter  Fji  (Table  10).  This  parameter  can  be  determined  for 
other  fluids  also;  in  this  project  the  procedure  has  not  been  carried  out  for  CO2.  For  the  calculations 
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Ffl-  1  has  been  used  due  to  similarities  of  evaporating  water  (compare  constants  used  in  VDI, 
1994). 


Tab 


e  10:  The  fluid-dependent  parameter  Ffi 


Huitl 

F, 

Water 

1.00 

R-ll 

1.30 

R-12 

1.50 

R-13B1 

1.31 

R-22 

2.20 

R-113 

1.30 

R-114 

1.24 

R-152a 

1.10 

n2 

4.70 

Neon 

3.50 

For  the  constants  Cm  there  are  two  columns,  one  for  convective  boiling  and  one  for  nucleate 
boiling  (Table  11).  The  greater  value  of  etc js, convective  and  cccj-5, nucleate  calculated  by  Eq  6  has  to  be 
used. 


Table  11:  Constants  Cm  for  the  Kandlikar  equation 


Constant 

convective 

boiling 

nucleate 

boiling 

Remark 

Cj 

1.1360 

0.66B3 

c2 

-0.9 

-0.2 

c3 

667.2 

1058.0 

c4 

0.7 

0.7 

c5 

0 

0 

vertical  tubes  for  any  Fr' ;  for  horizontal  tubes  with 
Fr'  >  0.04 

c5 

0.3 

0.3 

for  horizontal  tubes  with  Fr"  <  0.04 

One  has  to  consider  that  this  model  gave  only  acceptable  results  for  low  mass  fluxes. 


A.2.2  Pressure  Drop 

Thome  (1997)  recommends  the  Friedel  (1979)  correlation  for  the  frictional  pressure  gradient: 


f dp ^ 

Tl 


Jf  L 


2  -Z-m'Y 


n 


Eq  11 


where  $  is  a  two-phase  multiplier,  and  is  a  friction  factor  (see  Eq  19).  The  two-phase 
multiplier  is  defined  as 


4>'/=G  + 


3.2AF-H 


Fr 


0.045 


We 


0.035 


Eq  12 
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with  the  Froude  Number  (Frh)  and  the  Weber  Number  (Weh)  for  a 
without  slip  (Eq  16  and  Eq  17),  and  the  factors 

(fictitious)  homogeneous  flow 

G  =  (l-x)2  +x2  vAFl  , 
v  -%o 

Eq  13 

F  =  x078  \ l-*r',and 

Eq  14 

HJp'V  |YYi?  L  <Y7 

Ip  J  UJ  l 

Eq  15 

where  E,'0  and  E,"  are  given  by  Eq  19. 

The  Froude  Number  (Frh)  and  the  Weber  Number  (We/,)  for  a  homogeneous  flow  are  defined  as 

Fr= 

*  gd'Pk 

Eq  16 

We „  =  ’h‘-d 

O  -Pa 

Eq  17 

where  o  is  the  surface  tension  and  p /,  is  the  homogeneous  density: 


P* 


1— x 


The  Fanning-friction  factor  can  be  calculated  according  to: 


_  0.079 
~  i[R^ 


Eq  18 


Eq  19 


In  Eq  11  and  Eq  13,  the  subscript  ‘0’  refers  to  a  modified  Reynolds  Number  and  friction  factor, 
respectively,  which  have  to  be  calculated  with  the  properties  of  the  liquid  O  or  the  vapour  (") 
phase  and  with  the  total  mass  flux  (m ),  i.e.  assuming  a  single-phase  flow  at  the  same  mass  flux. 


For  a  two-phase  flow,  the  pressure  change  due  to  acceleration  becomes  according  to  the  theorem 
of  momentum: 


A pa  =  —  \m' -w  +  M* •W'X*  ~  m2  • 
A 


x2  ,  (*-*)2 
P'*e  P'M) 


Eq  20 


where  8  is  the  void  fraction.  Thome  (1997)  uses  the  approach  of  Premoli  (1970): 
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i+[s±2-4 

x  p 


Eq  21 


where  S  is  the  slip  ratio 


S  =  l+Er  77 — T~y'Ei 
i +y-E2 


Eq  22 


and  the  volumetric  vapour  fraction: 


The  Reynolds  Number  for  an  assumed  liquid  flow  is  defined  as 


where  m  is  the  total  mass  flux.  The  Weber  Number  is  defined  as 


Eq  23 


E,=  1.578  ■Re0~019-  ^ 


Eq  24 


E2  =  0.0273  ■  We  ■  Re0~°51  ■  ^ 


Eq  25 


P  P  -x 

•  x  +  p*  •(! — x) 


Eq  26 


,  m-d 
Reo=~ — 


Eq  27 


We  = 


m2  • d 


Eq  28 
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